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Nanoscale  magnetic  materials  exhibit  a wide  range  of  fascinating  phenomena 
which  have  applications  in  recent  technology.  Some  of  them  are  giant  magne- 
toresistance in  magnetic  multilayers,  Coulomb  blockade  in  tunnel  junctions,  and 
magnetoelectric  effect  in  spin-ordered  materials.  Understanding  the  physics  behind 
all  these  different  mechanisms  is  a challenge  to  any  theorist. 

Magnetic  multilayers  are  thin  alternating  layers  of  ferromagnets  and  paramag- 
nets.  In  these  multilayers  a large  change  in  resistivity  from  antiparallel  to  parallel 
alignment  of  magnetizations  hcis  been  observed.  This  new  phenomenon  is  called 
the  giant  magnetoresistance  and  is  being  applied  in  the  recording  industry.  The 
physics  of  giant  magnetoresistance  in  magnetic  multilayers  has  been  known  for 
several  years.  With  improved  fabrication  techniques  it  is  now  possible  to  fabricate 
nanowires  which  are  laterally  confined  magnetic  multilayers.  It  is  then  natural 
to  investigate  how  the  giant  magnetoresistance  changes  with  lateral  confinement. 
Theoretical  investigation  on  these  laterally  confined  magnetic  multilayers  revealed 
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that  in  general  the  giant  magnetoresistance  decreases  with  confinement.  This  re- 
duction can  be  understood  in  terms  of  decrease  in  the  effective  mean  free  path.  Also 
I find  that  there  are  special  parameter  regimes  where  the  giant  magnetoresistance 
increases  with  lateral  confinement. 

Transport  properties  in  tunnel  junctions  with  magnetic  metals  in  the  Coulomb 
blockade  regime  are  being  studied.  In  this  regime  the  current-voltage  curves  take 
the  form  of  a “staircase”  structure.  To  understand  the  experimental  results  we 
study  the  magnetoresistance  in  these  tunnel  junctions  and  in  the  Coulomb  blockade 
regime.  My  study  shows  that  the  magnetoresistance  curves  exhibit  spikes  as  a 
function  of  voltage.  The  spikes  occur  at  those  places  where  the  current  increases 
by  a step.  The  spikes  finally  disappear  at  large  voltages. 

Change  of  magnetization  (or  polarization)  by  applying  an  electric  field  (or  mag- 
netic field)  in  some  spin-ordered  materials  is  known  as  the  magnetoelectric  effect 
and  has  been  observed  experimentally  in  antiferromagnetic  insulators.  For  thin 
magnetic  films  this  effect  has  not  been  measured.  With  a simple  atomic  calculation 
I show  that  the  electric  field  can  indeed  change  the  orientation  of  magnetization 
by  changing  the  magnetic  surface  anisotropy  of  the  material. 
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CHAPTER  1 
INTRODUCTION 


Magnetic  nanostructures  are  electrical  devices  with  smallest  dimension  of  or- 
der of  a nanometer  (0.000000001  meters)  [1,  2].  Recent  advances  in  fabrication 
techniques  such  as  optical,  x-ray  and  electron-beam  lithography,  molecular-beam 
epitaxy  (MBE),  sputtering,  and  others  have  made  it  possible  to  fabricate  such 
nanostructures.  Magnetic  nanostructures  now  play  an  important  role  in  the  com- 
mercial market  as  well  as  in  scientific  research  activities.  Areas  of  application  in- 
clude the  automotive  and  aerospace  industries,  medicine,  surveillance,  information 
storage  systems,  communication  and  manufacturing  control.  Study  of  electrical 
transport  properties  thus  has  become  essential.  In  my  thesis  I discuss  the  trans- 
port properties  of  three  different  types  of  nanostructures:  laterally  confined  mag- 
netic multilayers,  double-tunnel-junction  with  magnetic  metals,  and  thin  magnetic 
films.  Before  detailing  the  specific  models  used  to  study  the  transport  properties 
of  the  above  mentioned  nanostructures  it  is  useful  to  review  the  general  methods, 
the  experimental  results,  and  their  applications. 

1.1  Giant  Magnetoresistance  (GMR)  in  Magnetic  Multilayers 

Magnetic  multilayers  are  thin  alternating  layers  of  ferromagnets  (FM)  and  para- 
magnets  (PM)  of  thicknesses  a few  angstroms  (15-30  A).  Fig.  1.1  is  a prototype 
of  a magnetic  multilayer.  In  these  layers  a large  decrease  in  the  resistivity  from 


1 


2 


antiparallel  to  parallel  alignment  of  the  film  magnetizations  has  been  observed 
experimentally.  This  new  phenomenon  is  known  as  the  giant  magnetoresistance 


Figure  1.1:  Schematic  diagram  of  a 3-layer  magnetic  structure  (FM/PM/FM).  A 
paramagnetic  spacer  layer  (PM)  is  sandwiched  between  two  ferromagnetic  (FM) 
layers.  The  direction  of  the  current  is  shown  with  the  arrows.  The  current  can  be 
in-plane  (CIP)  or  perpendicular  to  the  plane  (CPP). 

(GMR).  Although  there  were  earlier  reports  of  unusual  magnetoresistive  eflFects  in 
layered  structures  [3,  4,  5]  it  was  only  around  1988  that  GMR  effects  were  first 
clearly  observed  and  characterized  in  antiferromagnetically  coupled  Fe/Cr  systems. 
In  Fig.  1.2  we  show  the  first  observation  of  GMR  in  Fe/Cr  multilayers  by  Baibich 
et  al.  [6,  7]  and  one  observes  that  the  resistivity  drops  dramatically  when  an  ap- 
plied field  aligns  the  magnetic  moments  of  successive  layers.  Since  then,  the  quest 
for  larger  GMR  has  continued  with  different  layered  structures  and  films  grown 
either  by  molecular  beam  epitaxy  (MBE)  or  by  sputtering  techniques.  All  these 
efforts  lead  to  CIP  (current-in-plane)  - GMR  ranging  from  4%  - 150%  [8]  but 
also  have  provided  us  with  valuable  information.  Systematic  studies  of  GMR  as  a 
function  of  the  thickness  of  the  spacer  layer  have  shown  oscillations  in  the  mag- 
netoresistance. Fig.  1.3  shows  one  of  the  result  for  Co/Cu  multilayers  obtained 
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Figure  1.2:  Magnetoresistance  of  three  Fe/Cr  superlattices  at  4.2  K for  different 
Cr  thicknesses  as  measured  by  Baibich  et  al.  [6]  The  current  and  the  applied  field 
are  along  the  same  [110]  axis  in  the  plane  of  the  layers.  H,  is  the  saturation  field 
required  to  align  the  magnetizations  of  the  Fe  layers  in  parallel. 
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by  Parkin  [9,  10,  11].  A large  GMR  is  present  when  the  layers  are  antiferro- 


Figure  1.3:  Magnetoresistance  (MR)  for  Co  (10  A)/Cu  (tcu)  measured  by  Parkin  et 
al.  [9]  Magnetoresistance  oscillates  as  a function  of  the  spacer  layer  thickness  which 
indicates  that  there  is  an  interlayer  exchange  coupling  between  the  ferromagnetic 
layers. 

magnetically  coupled  and  a very  small  GMR  otherwise.  The  oscillations  indicate 
that  there  is  an  interlayer  exchange  coupling  of  Ruderman-Kittel-Kasuya-Yosida 
(RKKY)  type  between  the  magnetic  layers  [12].  The  GMR  is  generally  due  to  the 
interplay  between  spin  dependent  scattering  in  successive  magnetic  layers.  The 
conditions  for  the  interplay  is  that  the  distance  between  the  layers  is  relatively 
small  in  comparison  with  the  mean  free  path  of  the  electron.  The  variation  of  the 
GMR  with  the  thickness  of  the  magnetic  layers  also  contains  potential  information 
about  the  ratio  of  bulk  and  interfacial  spin  dependent  scattering.  In  this  context 
it  is  useful  to  mention  that  bulk  scattering  means  that  the  mean  free  paths  of  the 
electrons  are  spin-dependent.  On  the  other  hand  interface  scattering  arises  due  to 
the  difference  between  transmission  and  reflection  coefficients  for  the  spin-up  and 
spin-down  electrons  at  the  interface. 

Influence  of  interface  structure  plays  an  important  role  on  the  GMR  and  also 
have  been  studied  by  various  groups  [8,  13,  14].  It  has  been  found  that  the  spin- 
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dependent  scattering  by  the  interfaces  is  enhanced  when  the  surfaces  are  “rougher” 
which  thus  increases  the  GMR. 

GMR  also  depends  on  the  direction  of  the  current.  In  most  experiments  the 
current  flows  in  the  plane  of  the  layers,  this  is  the  CIP  (current-in-plane)  geometry. 
However,  magnetoresistance  measurements  also  have  been  done  for  the  case  where 
the  current  is  perpendicular  to  the  plane  - this  is  known  as  the  GPP  geometry  (see 
Fig.  1.1).  Vertical  transport  (GPP)  is  the  common  mode  for  using  semiconducting 
superlattices,  but  for  a given  sample  of  1 //m  thick  and  area  1 mm^  in  the  plane  of 
the  layer  the  resistance  is  one  million  times  smaller  for  GPP  than  CIP.  Thus,  it  is 
harder  to  measure  the  GPP  resistance  by  conventional  means.  However,  CPP-MR 
ratio  has  been  measured  in  sputtered  Co/Cu  and  Ag/Co  superlattices  by  a group 
in  Michigan  State  University  [15,  16]  and  a value  of  170%  has  been  obtained.  The 
GMR  in  GPP  geometry  exhibits  similar  oscillations  as  a function  of  the  thickness 
of  the  non-magnetic  spacer  layers,  and  decreases  less  rapidly  at  large  thicknesses. 
Fig.  1.4  shows  the  experimental  result  [15,  16]. 

Temperature  dependence  of  the  GMR  is  of  critical  importance  to  its  applica- 
tions. It  has  been  seen  that  the  MR  ratio  drops  with  increase  in  temperature  [13]. 
Fig.  1.5  illustrates  this  feature.  Several  aspects  contribute  to  the  decrease  in  GMR. 
At  higher  temperature  inelastic  scattering  processes  (with  phonons,  magnons,  etc.) 
come  into  play  which  reduces  the  mean  free  path  which  then  decreases  the  mag- 
netoresistance. 

Since  the  discovery  of  giant  magnetoresistance  in  magnetic  multilayers  many 
different  theoretical  models,  both  classical  and  quantum,  have  been  proposed  to 
understand  the  various  features  of  the  GMR  seen  in  the  experiment  [17,  18,  19,  20]. 
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Figure  1.4:  The  magnetoresistance  of  Co/Cu  layers  at  T =4.2  K as  a function  of 
the  thickness  of  the  spacer  layer  (Cu)  for  currents  perpendicular  to  the  plane  of  the 
layers  (CPP)  as  measured  by  the  Michigan  State  University  group.  The  change 
in  resistance  is  measured  from  the  zero  field  value  to  saturation  (Hq),  and  from 
its  peak  value  after  recycling  through  the  hysteresis  curve  (Hp).  Comparing  the 
reproducible  peak  values  of  the  CIP-MR  with  CPP-MR  shows  that  CPP-MR  is 
much  larger  than  the  CIP-MR. 
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Figure  1.5:  Magnetoresistance  as  a function  of  the  magnetic  field  at  different  tem- 
peratures as  measured  by  Bader’s  group  [8].  A maximum  CIP-MR  of  150%  is 
observed  at  4.2  K (28%  at  room  temperature)  for  Fe/Cr  superlattices. 
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Fuchs-Sondheimer’s  semiclaissical  theory  for  the  conduction  in  thin  films  has  been 
known  since  1950s  [21,  22].  It  was  later  extended  by  Camley  and  Barnas  for  the  CIP 
geometry  to  magnetic  multilayers  [23].  Their  model  has  been  used  extensively  to 
numerically  calculate  the  magnetoresistance  in  multilayers  [24,  25,  26]  . In  the  limit 
where  the  thickness  of  the  magnetic  multilayers  is  much  larger  than  the  mean  free 
path  of  the  conduction  electrons  Barthelemy  and  Fert  obtained  simple  analytical 
expressions  for  the  GMR  [27].  In  the  limit  of  very  long  mean  free  path  analytical 
results  have  been  obtained  by  Edwards  et  al.  [28]  Various  modifications  of  the 
model  of  Camley  and  Barnas  have  been  proposed.  For  example,  to  describe  the 
scattering  at  the  interface  Johnson  and  Camley  [29]  introduced  a thin  interfacial 
layer  with  a shorter  mean  free  path.  Hood  and  Falicov  took  into  account  the 
potential  scattering  by  the  interfaces  [30].  For  CPP  geometry  theoretical  models 
based  on  Boltzmann  equation  have  also  been  worked  out  [31,  32]. 

In  the  regime  when  the  film  thickness  becomes  smaller  than  the  electron  mean 
free  path,  the  semiclassical  theory  of  Fuchs  and  Sondheimer  predicts  that  the 
resistivity  varies  as  the  inverse  of  the  logarithm  of  the  mean  free  path.  In  the 
limit  of  infinite  mean  free  path,  the  resistivity  becomes  zero  which  implies  that  in 
the  absence  of  bulk  scattering  there  is  no  contribution  due  to  surface  scattering. 
One  would  rather  expect  a finite  surface  scattering  and  therefore  a finite  residual 
resistivity  in  this  limit.  To  correct  this  non-physical  result,  the  quantum  theory 
of  transport  in  thin  metallic  films  in  the  presence  of  disorder  was  first  developed 
by  Tesanovic  et  al.  [33]  This  can  also  be  understood  by  considering  the  Landauer 
conductance  formula  [34]  in  the  absence  of  disorder  and  for  temperatures  much 
less  than  the  Fermi  energy  {ksT  <<  Ep).  The  conductance  {G)  is  expressed  in 
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terms  of  the  number  of  channels  (N).  For  a 3D  film  of  length  L with  finite  cross 
section  S the  number  of  channels  N — Skp/4n  where  kp  is  the  Fermi  wave  vector. 
The  end  result  for  the  conductance  is 


aS 

~T 


(1.1) 


where  cr  is  the  conductivity.  The  conductivity  is  now  obtained  from  the  conduc- 
tance as 

This  result  shows  that  the  conductivity  is  always  finite  for  a sample  of  finite  size. 
Different  quantum  models  of  GMR  based  on  Kubo  formalism  [17,  35,  36,  37,  38, 
39,  40,  41],  Landauer-Buttiker  formalism  [42],  quantum  Boltzmann  equation  [43], 
ab-initio  calculations  using  realistic  superlattice  band  structures  [44,  45,  46,  47] 
etc.  have  been  developed. 

The  basic  qualitative  features  of  the  giant  magnetoresistance  can  be  understood 
classically  by  a simple  model  [19].  In  this  model  one  assumes  that  the  GMR  is 
entirely  due  to  the  bulk  scattering,  the  mean  free  path  of  the  electrons  within  the 
layers  (100  A is  a typical  value)  is  much  larger  than  the  thickness  of  the  layers, 
and  there  is  no  spin-flip  scattering.  Fig.  1.6  schematically  shows  the  mechanism 
of  the  magnetoresistance.  For  high  magnetic  fields  the  magnetizations  are  parallel. 
One  can  then  imagine  this  as  an  electrical  circuit  with  two  resistivity  channels 
(for  majority  and  minority  electrons)  in  parallel  (Fig.  1.6a).  Depending  on  the 
spin  directions  — 1/2  or  = —1/2  the  scattering  probabilities  and  therefore 
the  resistivities  are  different.  The  majority  electrons  with  = 1/2  form  a low 
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(a)  (b) 


Figure  1.6;  Schematic  diagram  of  conduction  in  multilayer  magnetic  films  showing 
how  spin  scattering  produces  a different  resistance  for  parallel  (case  a)  and  anti- 
parallel (case  b)  alignments  of  magnetizations  [19].  When  the  magnetizations  are 
parallel  (case  a)  it  is  relatively  easy  for  spin-up  (S^)  electrons  (majorities)  to 
move  freely  within  the  layers  (less  scattering)  where  as  spin-down  (S~)  electrons 
(minorities)  scatter  the  most.  Equivalent  resistor  array  for  parallel  alignments 
is  shown  below.  R and  r are  the  resistances  for  the  spin-down  and  the  spin-up 
electrons.  For  case  (a)  r is  smaller  than  R thus  the  effective  resistance  is  lower.  For 
case  (b)  magnetizations  are  anti-parallel  - no  electron  can  traverse  two  magnetic 
layers  without  becoming  the  unfavored  highly  scattered  species.  Thus  in  this  case 
the  high  and  the  low  resistances  are  in  series  and  the  effective  resistance  is  larger 
than  the  eflFective  resistance  in  case  (a). 
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resistivity  channel  throughout  the  sample. 

_ 

(p'+p^y 


(1.3) 


Fig.  1.6b  represents  the  opposite  case  with  an  antiparallel  alignments  of  the  mag- 
netizations in  successive  layers.  In  this  case  the  low  resistivity  electrons  (majority 
electrons)  in  a layer  becomes  the  high  resistivity  electrons  (minority  electrons) 
in  the  next.  Each  channel  hais  the  average  resistivity  {p^  + P'*')/2  and  the  final 
resistivity  is  high 


Paf 


(P^  + P^) 
4 


(1.4) 


The  magnetoresistance  in  the  infinite  mean  free  path  limit  is  then: 


GMR  = ~ 

Paf 


pi  + ptj 


(1.5) 


The  above  result  shows  that  in  the  absence  of  bulk  scattering,  i.e.  there 

is  no  magnetoresistance.  But  this  is  not  true  for  a real  system.  GMR  is  due 
to  both  bulk  and  interface  scattering,  so  even  in  the  absence  of  bulk  scattering 
there  is  GMR  due  to  the  interface  scattering  alone.  The  above  simple  model  can 
be  generalized  for  a FM/PM/FM  multilayer  with  different  thicknesses.  Let  tpM 
and  tpM  are  the  thicknesses  of  the  ferromagnet  and  the  paramagnet,  respectively. 
Then  the  total  resistivity  for  the  majority  electrons  (with  Sz  = 1/2)  is  obtained  by 
adding  the  resistivities  for  the  individual  layers  in  series.  In  terms  of  the  resistivity 


12 


(p’’^)  of  the  majority  electrons  we  have: 


ft  2tpMP^  + ^pmP 

p' ' - 

2tpM  + ^PM 


PM 


(1.6) 


Similarly,  the  resistivity  (p-*^)  for  the  minority  electrons  (with  = —1/2)  is  given 
by: 


2tFMP'*'  + tpuP 
2tpM  + ^PM 


PM 


(1.7) 


The  total  resistivity  for  the  parallel  alignments  of  the  magnetizations  is  then  ob- 
tained by  adding  the  resistivities  in  Eqs.  (1.6)  and  (1.7)  in  parallel: 


pttpU 

f’'  ~ ptt  + 


(1.8) 


Proceeding  in  the  same  way  the  total  resistivity  for  the  anti-parallel  alignments  of 
the  magnetizations  (in  this  case  is 


PAF 


ptj.  -\-  piA 

T” 

^fmP^  + ^pmP^^  + ^fmP'*' 
2(2tpM  + ^pm) 


(1.9) 


Using  Eq.  (1.5)  the  magnetoresistance  for  bulk  scattering  is  now  given  by: 


GMR(Bulk) 


t-eu{p^  - P^) 
^fm(p^  + P^)  + ^pmP^^ 


(1.10) 


In  the  limit  of  tpM  ->  0 we  recover  Eq.  (1.5)  from  Eq.  (1.10). 
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These  new  phenomenon  in  magnetism  has  been  combined  with  electronics  to 
build  new  devices  [48].  For  example,  magnetic  sensors  have  wide  applications  in 
our  present  day  technology  [49].  Applications  include  magneto-optic  sensors,  anti- 
theft devices.  Hall  effect  sensors,  different  kinds  of  magnetometers  such  as  fiber- 
optic magnetometers,  SQUID  magnetometers,  search-coil  magnetometers,  flux- 
gate  magnetometers,  and  many  others.  The  motivation  driving  the  research  in 
GMR  both  theoretically  and  experimentally  is  its  application  in  the  hard  disk  drive 
industry  [50].  High  sensitive  read  heads  made  of  giant  magnetoresistive  materials 
have  been  introduced  recently.  The  GMR  head  employs  a sensor  whose  resistance 
changes  in  the  presence  of  a magnetic  field.  Its  performance  gain  has  enhanced  the 
density  of  storage  by  more  than  50%.  GMR  also  increase  the  MR  strip  sensitivity 
dR/R  to  provide  higher  signal  amplitudes  [51,  52].  A prototypical  MR  head  is 
shown  in  Fig.  1.7.  The  head  senses  the  magnetic  field  from  the  recorded  bit.  With 


Figure  1.7:  Schematic  diagram  of  a magnetoresistive  “read”  head  [50].  With  a 
GMR  material  in  the  “read  head” , the  resistance  changes  with  a small  change  in 
the  magnetic  field. 


a GMR  material  in  the  sensor  the  resistance  of  the  material  changes  with  change  in 
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the  magnetic  field.  For  better  performance  one  needs  larger  sensitivity  i.e.  dR/dH 
to  be  large  which  means  a large  change  in  the  resistance  with  a small  change  in 
the  magnetic  field. 

Another  application  of  GMR  materials  is  in  the  magnetic  random  access  mem- 
ory (MRAM)  [53,  54].  In  MRAM  technology  two  particular  properties  of  mag- 
netic materials  are  used:  uniaxial  anisotropy  and  magnetoresistance.  Uniaxial 
anisotropy  aligns  the  magnetization  along  the  easy  axis.  The  magnetization  will 
thus  point  in  one  of  the  two  antiparallel  directions  which  are  used  to  store  “1”  or 
a “0” . The  magnetoresistance  property  of  the  material  can  then  be  used  to  distin- 
guish between  a “1”  or  a “0”  in  a MRAM  cell.  The  basic  MRAM  cell  configuration 
is  a magnetic  sensing  line  and  an  orthogonal  and  insulated  word  line  which  is  shown 
in  Fig.  1.8.  The  magnetic  sense  line  consists  of  thin  magnetic  films  sandwiched 


Figure  1.8:  A GMR  based  MRAM  cell:  (a)  cross-section  along  the  sense  line  and 
(b)  the  direction  of  the  magnetic  field  across  the  sense  line  due  to  a current  in  the 
line  which  goes  into  the  paper  [53]. 


between  a nonmagnetic  layer.  The  magnetic  thin  films  are  made  of  giant  mag- 
netoresistive materials.  Due  to  interlayer  exchange  coupling  the  magnetic  fields 
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align  across  the  sense  line  in  opposite  directions  in  the  top  and  the  bottom  film  (as 
shown  by  the  arrows  in  Fig.  1.8b)  with  the  application  of  a current  (which  goes 
into  the  page).  About  the  sense  current  the  storage  states  are  either  in  clockwise 
or  in  anticlockwise  orientation  of  the  magnetizations  in  the  sandwich  layers.  A one 
or  zero  is  written  into  a two-dimensional  array  of  these  cells  by  a word  current  in 
combination  with  a sense  current  of  appropriate  polarity.  MRAM  chips  with  up 
to  16  kilobits  of  storage  have  been  demonstrated.  Designs  of  chips  with  up  to  a 
million  bits  are  now  in  progress. 

With  such  a tremendous  progress  in  nanofabrication  techniques  it  is  now  pos- 
sible to  laterally  confine  the  magnetic  multilayers  [55,  56,  57].  Multilayered  wire 


Figure  1.9;  Schematic  diagram  of  a laterally  confined  magnetic  multilayer.  PM  is 
the  paramagnetic  layer  sandwiched  between  two  ferromagnetic  (FM)  layers.  The 
direction  of  the  electric  field  (E)  is  in  the  plane  of  the  layers  (CIP). 

arrays  have  been  fabricated  recently  by  high-resolution  electron-beam  lithogra- 
phy and  lift-off  techniques  to  study  finite  size  effects  and  giant  magnetoresis- 
tance [58,  59].  A prototype  of  a laterally  confined  magnetic  multilayer  is  shown 
in  Fig.  1.9.  Laterally  confined  magnetic  multilayers  hold  the  potential  to  be  used 


Width 
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in  future  technologies  [60].  It  is  then  natural  to  study  these  laterally  confined 
magnetic  multilayers.  Understanding  the  effect  of  lateral  confinement  on  GMR 
theoretically  is  the  primary  focus  of  my  investigation  [61]. 

1.2  Junction  Magnetoresistance  (JMR)  and  Coulomb  Blockade 

in  Tunnel  Junctions 

Electron  tunneling  is  a quantum  phenomenon  in  which  the  electrons  travel  from 
one  electrode  to  another  through  a thin  oxide  layer  which  acts  as  an  insulating  bar- 
rier. This  three-layer  system  - electrode,  barrier,  and  counter-electrode  is  referred 
to  as  tunnel  junction.  The  quantity  usually  measured  in  a tunneling  experiment 
is  the  current  or  its  derivative  as  a function  of  applied  voltage  [62]. 

Some  of  the  earliest  experiments  were  on  tunneling  between  a superconductor 
and  a ferromagnet  [63].  In  these  experiments,  the  difference  between  the  spin- 
up  and  spin-down  density  of  states  of  a superconductor  in  a magnetic  field  was 
observed  [64] . These  experiments  also  provided  a quantitative  measure  of  the  spin 
polarization  of  the  current  from  the  ferromagnetic  metal. 

Later  experiments  were  done  on  tunneling  between  two  ferromagnetic  met- 
als [65].  Because  of  the  spin  dependence  of  the  density  of  states,  the  conductance 
changed  when  the  relative  orientation  of  the  two  magnetizations  changed.  One 
of  the  measurable  quantities  of  practical  importance  is  the  junction  magnetoresis- 
tance, which  is  defined  as  the  ratio  of  the  change  in  the  conductance  from  parallel 
to  antiparallel  alignment  divided  by  the  parallel  conductance.  Large  magnetore- 
sistance response  in  low  magnetic  fields  has  been  observed  in  a variety  of  tunnel 
junctions  by  several  groups  [66,  67,  68,  69,  70,  71,  72,  73,  74,  75,  76,  77].  Fig.  1.10 
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Figure  1.10:  Magnetoresistance  for  NiMnSb/Al203/NiFe  junction  at  room  tem- 
perature and  77  K [66].  The  left  scale  corresponds  to  the  room  temperature  (RT) 
curve,  the  right  to  the  77  K data.  The  broad  peak  at  77  K indicates  well-separated 
coercive  forces  of  the  ferromagnets  and  antiparallel  orientation  at  the  peak.  The 
sharp  peak  at  RT  indicates  that  the  coercivities  of  the  two  FM’s  are  close  to  each 
other  which  implies  that  complete  antiparallel  orientation  is  not  achieved  at  the 
peaks.  The  maximum  value  obtained  for  JMR  is  about  2.4%  at  RT  where  as 
highest  JMR  observed  at  77  K is  about  7.0%. 
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shows  the  JMR  as  a function  of  the  magnetic  field  measured  by  Tanaka  et  al.  [66] 
From  the  point  of  view  of  application  higher  value  of  JMR  as  well  as  low  junction 
resistance  of  the  order  of  1 kO//im^  is  required.  Very  recently  such  a value  for  the 
resistance  has  been  obtained  [78].  Using  half-metallic  compounds  as  electrodes  can 
also  increase  the  value  of  JMR.  At  the  Fermi  level  the  half-metallic  compounds 
have  only  one  spin-direction  which  is  reflected  as  a gap  in  the  density  of  states. 
Examples  of  such  compounds  are  Heusler  alloys  [79,  80,  81]. 

To  estimate  the  junction  magnetoresistance,  Julliere  in  1975  proposed  a sim- 
ple model  based  on  the  free  electron  picture  [65].  For  a ferromagnetic-insulator- 
ferromagnetic  (FM/I/FM)  tunnel  junction  the  conductance  depends  on  the  relative 
magnetizations  of  the  two  ferromagnetic  electrodes.  For  magnetizations  parallel  in 
both  the  electrodes  we  define  the  conductance  to  be  and  is  the  conduc- 
tance when  the  magnetization  of  the  electrodes  are  anti-parallel.  Let  oi  and  02  be 
the  fraction  of  the  number  of  spins  parallel  to  the  magnetization  and  N\  and  N2 
the  density  of  states  of  the  electrons  at  the  Fermi  level  in  the  two  ferromagnets. 
Using  Fermi’s  golden  rule  we  obtain  the  two  conductances  which  are 

G^  = |T|2iVxiV2  {aia2  + (1  - ai)(l  - 02)}  , (1-11) 

G^^  = \T\'^NrN2  {ai(l  - 02)  + (1  - ai)a2>  . (1.12) 

Here  the  tunneling  matrix  element  T,  is  taken  to  be  both  energy  and  momentum 
independent.  Using  the  definitions  of  the  spin  polarizations  Pi  and  P2  of  the  two 
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ferromagnets 


Pi  = 2ax  - 1,  (1.13) 

P2  = 202-1.  (1.14) 


Eqs.  (1.11)  and  (1.12)  can  be  rewritten  in  terms  of  Pi  and  P2.  The  junction  mag- 
netoresistance (JMR)  is  defined  as  the  ratio  of  the  change  in  the  conductance  from 
parallel,  to  antiparallel  alignments,  divided  by  the  parallel  conductance 

^AP 

JMR  = 1--^.  (1.15) 

The  above  expression  for  JMR  can  also  be  rewritten  in  terms  of  the  polarizations 


JMR  = 


2P1P2 

(1  + AA)' 


(1.16) 


For  a Fe/Ge/Co  junction  using  the  polarization  values  (Pi=44%,  P2=34%)  as 
measured  by  Tedrow  and  Meservey  on  Al-Al203-ferro  junctions  Julliere  estimated 
the  value  of  JMR  from  Eq.  (1.16)  to  be  26%.  The  maximum  measured  experi- 
mental value  for  this  Fe/Ge/Co  tunnel  junction  is  about  14%.  The  discrepancy 
between  the  experimental  value  and  the  theoretically  predicted  value  may  be  due 
to  spin-flip  processes,  magnetic  impurities,  and  other  mechanisms. 

Besides  spin,  electrons  carry  charge  and  at  low  temperatures  in  tunnel  junctions 
with  small  capacitance,  the  discreteness  of  electron  charge  manifests  itself  through 
the  physics  of  the  Coulomb  blockade  [82,  83,  84,  85,  86,  87].  The  transfer  of 
an  electron  by  tunneling  between  two  initially  neutral  regions  of  capacitance  C 
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increaises  the  electrostatic  energy  of  the  system  by  an  amount  e^/2C.  At  low 
temperatures  and  small  voltages,  the  tunneling  current  is  suppressed  because  of 
the  charging  energy.  This  energy  barrier  is  called  the  Coulomb  blockade.  The  effect 
of  this  in  a two-junction  system  (a  prototype  of  which  is  shown  in  Fig.  1.11)  is  the 
incremental  increase  in  the  current  at  voltages  where  it  is  energetically  favorable 
for  an  electron  to  tunnel  to  the  center  island.  The  occurrence  of  these  current  steps 
in  the  current-voltage  characteristics  is  known  as  the  “Coulomb  staircase.”  One 


Figure  1.11;  The  basic  structure  of  a double  tunnel  junction.  The  source,  island, 
and  the  drain  are  the  electrodes  separated  by  insulating  layers. 

of  the  recent  experiments  on  the  Coulomb  blockade  phenomenon  in  nanoparticles 
is  by  Ralph  et  al.  [88]  Using  electron-beam  lithography  and  reactive-ion  etching 
techniques  his  group  fabricated  a bowl-shaped  hole  in  a suspended  silicon  nitride 
membrane  (See  Fig.  1.12a).  The  orifice  has  a size  between  5 and  10  nm  in  diameter. 
A single  particle  is  then  formed  by  deposition  techniques  under  the  nm-scale  tunnel 
junction  to  contact  the  top  and  the  bottom  Al  electrodes.  Electrons  tunnel  from 
the  top  to  the  bottom  electrode  through  the  single  nano-particle.  The  gate  voltage 
tunes  the  number  of  electrons  in  the  particle.  In  this  device  geometry  the  current 
plotted  as  a function  of  the  voltage  clearly  shows  the  Coulomb  staircase  structure 
(Fig.  1.12). 
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Figure  1.12:  (a)  Schematic  cross  section  of  the  device  geometry.  Electrons  tunnel 
between  the  top  and  the  bottom  electrodes  through  a single  nanoparticle.  The 
gate  voltage  Vg  controls  the  flow  of  the  electrons.  Figs,  (b)  - (d)  display  the 
Coulomb-staircase  structure  for  different  samples,  at  equally  spaced  values  of  gate 
voltage  [88]. 
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There  are  number  of  new  experiments  which  exhibit  the  Coulomb  blockade 
phenomenon:  granular  thin  films  [89,  90],  lithographically  patterned  tunnel  junc- 
tions [88,  91,  92,  93,  94],  narrow  disordered  quantum  wires  [95,  96],  scanning 
tunneling  microscopy  of  small  metal  droplets  [97],  metal-insulator-metal  tunnel 
junctions  with  small  metal  particles  embedded  in  the  insulator  [98,  99,  100,  101]. 

Recently  there  have  been  experiments  where  spin-polarized  tunneling  and  the 
Coulomb  blockade  have  been  combined  [102,  103,  104,  105,  106].  In  discontinuous 
metal/insulator  multilayers  consisting  of  closely  spaced  ferromagnetic  nanoparti- 
cles in  an  insulating  matrix  a Coulomb  blockade  is  seen  (Fig.  1.13)  in  the  current- 
voltage  characteristic  at  low  temperatures. 


Figure  1.13:  I(V)  at  different  temperatures  for  the  annealed  multilayer  [CoFe(15 
A)HfO2(40  A)]2o  [102].  At  low  temperature  we  see  a Coulomb  blockade.  In  the 
inset  the  differential  conductance  (dl/dV)  at  10  K is  plotted  to  show  the  magnitude 
of  the  Coulomb  blockade  (0.3  V). 
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To  explain  the  recent  experimental  results  theories  have  been  proposed  [107, 
108,  109].  The  role  of  electronic  structure  in  magnetic  tunneling  has  also  been 
studied  [110,  111,  112], 

Spin-transistors  and  spin-polarized  field-effect-transistor  (known  as  spin  FET) 
are  based  on  the  application  of  spin-polarized  transport  [113].  In  the  near  future 
it  may  be  possible  to  make  a device  which  combines  both  of  these  effects  (spin- 
polarized  tunneling  and  Coulomb  blockade)  [114].  The  purpose  of  my  work  is  to 
theoretically  study  the  different  transport  properties  of  a double-tunnel-junction 
with  magnetic  metals  at  low  temperatures  where  the  phenomenon  of  Coulomb 
blockade  ha.s  been  seen  [108]. 

1.3  Magnetization  Reversal  in  Ultra-Thin  Magnetic  Films  and 

in  Ferromagnetic  Grains 

One  of  the  present  day  research  areas  in  magnetics  is  focussed  on  the  phe- 
nomenon of  magnetization  reversal  in  ultra-thin  magnetic  films  and  in  single  do- 
main magnetic  nano-particles  [115,  116,  117,  118,  119].  Magnetization  reversal  in 
nanoscale  ferromagnets  has  potential  application  in  the  next-generation  recording 
media  such  as  magnetic  tapes  and  disks.  During  recording  different  regions  of 
the  recording  medium  are  briefly  exposed  to  strong  magnetic  fields  so  that  each 
grain  is  magnetized  in  a certain  direction.  For  better  storage  devices  the  grains 
should  retain  their  magnetization  for  long  periods  of  time.  The  grain  size  and  the 
lifetime  of  the  magnetization  opposed  to  the  applied  magnetic  field  are  therefore 
of  considerable  technological  interest. 
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Magnetization  reversal  by  droplet  formation  is  a possible  scenario  in  single- 
domain nanoscale  ferromagnets  [120,  121,  122,  123,  124,  125,  126,  127]  but  the 
standard  theory  is  due  to  Neel  [128]  and  Brown  [129,  130].  Before  we  discuss  Neel- 
Brown’s  theory  [131]  let  us  illustrate  the  phenomenon  of  magnetization  reversal 
with  a simple  example  from  statistical  mechanics.  Consider  a spin-chain.  We  ap- 
ply a magnetic  field  large  enough  to  saturate  the  magnetization  of  this  chain  in  the 
field’s  direction.  We  then  slowly  reduce  the  field  to  zero  and  then  increase  slowly  in 
the  opposite  direction.  At  some  stage  during  this  process,  the  initial  state  will  be- 
come unstable,  and  some  change  will  start  to  take  place.  Finally  the  magnetization 
will  saturate  in  the  opposite  direction.  The  field  at  which  the  original  saturated 
state  becomes  unstable  and  any  change  in  magnetization  just  starts  is  called  the 
nucleation  field  (also  known  as  the  switching  field).  The  transition  state  between 
the  initial  and  the  final  stable  states  is  known  as  the  metastable  state.  Metastable 
phases  are  observed  in  a wide  variety  of  systems  in  condensed  matter  physics,  high 
energy  physics,  and  in  cosmology.  A few  examples  are  supercooled  liquids,  perma- 
nent magnets,  ferroelectrics,  in  certain  alloys,  in  the  electroweak  phase  transition, 
and  the  supercooled  quark-gluon  plasma  associated  with  the  QCD  confinement 
transition  [132].  A fully  satisfactory  description  of  the  metastable  phase  would  be 
extremely  helpful  to  understand  all  these  different  systems.  Some  of  the  pertinent 
questions  to  ask  are:  (1)  How  does  a metastable  state  decay  into  a more  stable 
state?  (2)  How  long  is  the  switching  time,  i.e.  the  time  to  go  from  one  stable  state 
to  another  stable  state?  (3)  What  are  the  factors  the  switching  time  depend  on 
and  how  can  one  control  the  switching  time? 
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Nanoscale  ferromagnetic  grains  have  been  studied  for  many  years,  but  recently 
with  magnetic  force  microscopy  (MFM)  and  Lorentz  microscopy  the  magnetic 
properties  of  individual  single-domain  particles  are  studied.  Chang  et  al.  [133]  used 
MFM  to  study  the  dependence  of  the  switching  field  with  the  diameter  of  isolated 
barium  ferrite  fine  particles.  The  resulting  effective  field  has  a peak  at  a grain 
diameter  of  about  55  nm.  MFM  measurements  in  individual  iron  particles  [134] 
showed  that  the  switching  fields  decrease  with  the  increase  of  the  particle  diameters 
from  20  to  70  nm.  Study  of  size  effects  in  single-domain  Ni  bars  [135]  also  shows 
peak  in  the  switching  field  at  a width  of  55  nm  which  then  decreases  with  increase 
in  the  bar-width.  Another  set  of  experiments  involved  the  measurement  of  the 
probability  P that  the  magnetization  in  7-Fe203  [119]  single-domain  particles  is 
switched  by  an  opposing  field.  The  experiments  show  that  the  mean  switching 
time  depends  strongly  on  the  strength  of  the  applied  field. 

To  understand  the  phenomenon  of  magnetization  reversal  we  need  to  study 
Brown’s  equation  of  micromagnetics  [130,  136].  We  consider  a ferromagnetic  body 
of  any  shape  in  which  the  magnetization  is  a function  of  space.  The  total  energy 
(e)  then  consists  of  four  terms:  the  exchange  energy  (cg),  the  anisotropy  energy 
(ca),  the  magnetostatic  energy  (cm),  and  the  Zeeman  energy  (e^): 


e — Cg  -f  Co  -h  Cm  + 6h-  (1-17) 

The  only  term  which  is  left  out  in  Eq.  (1.17)  is  magnetostriction.  When  a ferro- 
magnet  is  magnetized  it  shrinks  (or  expands)  in  the  direction  of  the  magnetization 
- this  is  know  as  magnetostriction.  In  some  cases  magnetostriction  can  be  absorbed 
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inside  the  anisotropy  term.  The  exchange  energy  stems  from  the  interaction  be- 
tween the  neighboring  spins.  Written  in  terms  of  the  magnetization  unit  vectors 
it  takes  the  following  form: 


where  C is  a constant  which  depends  on  the  coupling  constant  (assumed  to  be 
constant),  lattice  spacing,  and  size  of  the  spins. 

For  most  cases  the  Hamiltonian  for  a crystal  is  not  isotropic  which  means  that 
the  energy  levels  depend  on  the  direction  in  space  in  which  the  crystal  is  magne- 
tized. This  is  known  as  anisotropy  [137].  For  example  the  anisotropy  of  hexagonal 
crystals  depend  only  on  the  angle  (^)  between  the  c-axis  and  the  direction  of  mag- 
netization. The  anisotropy  energy  density  in  this  case  is 


where  Ki  and  K2  are  known  as  anisotropy  constants.  Depending  on  the  magnetic 
material  Ki  and  K2  can  have  both  signs.  K2  is  usually  an  order  of  magnitude 
smaller  than  Ki,  so  the  sign  of  Ki  determines  whether  the  magnetization  prefers 
in-plane  or  perpendicular  to  the  plane.  For  example  in  Fe  films,  Ki  is  positive  [138, 
139,  140,  141,  142],  so  the  magnetization  prefers  to  be  perpendicular  to  the  plane. 
On  the  other  hand,  Ki  is  negative  for  Ni  and  Co  films  - magnetization  thus  prefers 


(1.18) 


Wa  = Ki  sin^  0 + K2  sin^  9 

= Ki{l  - ml)  + K2{1  - mlf , 


(1.19) 
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to  be  in-plane  [140].  In  general  the  total  anisotropy  energy  can  be  written  as 


€a  = J dr  Wa  + j ds  w,.  (1.20) 

Here  Wa  and  Wg  are  the  volume  and  the  surface  anisotropy  densities,  respectively. 

The  origin  of  the  magnetostatic  term  is  the  coupling  of  the  magnetization  with 
the  internal  field  (Hint)  due  to  the  dipoles: 


Cm  = -^y^rM-Hint.  (1.21) 

The  remaining  Zeeman  term  is  due  to  the  response  of  the  external  magnetic  field 

(H.) 

€h  = - JdrM-H^.  (1.22) 

Brown’s  equation  can  be  obtained  by  minimizing  the  total  energy  in  Eq.  (1.17) 
and  the  result  is 

m X HeflF  = 0,  (1.23) 


with  the  constraint  |m|  = 1.  Hes  in  Eq.  (1.23)  is  defined  as 


Heff  = 


Mg 


V^m  + H - 


1 dWa 
Mg  dm 


(1.24) 


Mg  is  the  saturation  magnetization  in  the  above  equation.  Eq.  (1.23)  implies  that 
in  equilibrium  the  total  torque  is  zero  everywhere  and  the  magnetization  is  parallel 
to  an  effective  field.  Brown’s  equations  have  been  solved  for  different  geometries 
and  the  nucleation  fields  also  have  been  calculated  in  each  case.  Buckling,  curling. 
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fanning,  uniform  rotation  are  different  theoretical  relaxation  processes  for  different 
cases  [131],  Using  the  effective  field  we  can  also  study  the  details  of  magnetization 
configurations  and  magnetization  reversal  processes  in  magnetic  materials.  The 
dynamics  of  magnetization  orientation  follows  the  famous  Landau-Lifshitz-Gilbert 
equation: 

dM  A 

— = -7M  X H - — M X (M  X H),  (1.25) 

where  7 is  the  electron  gyromagnetic  ratio  and  A is  the  damping  constant.  Based 
on  Eq.  (1.25)  which  is  the  basic  equation  of  micromagnetic  theory  one  can  simulate 
complete  magnetization  processes  for  different  geometries  [143,  144,  145].  Micro- 
magnetic  modeling  enriches  our  understanding  about  the  formation  of  complex 
domains. 

In  Neel-Brown’s  theory  one  assumes  uniform  rotation  of  all  the  atomic  moments 
in  the  system.  This  is  to  avoid  the  energy  barrier  due  to  exchange  interactions 
between  atomic  moments  with  unlike  orientations.  The  remaining  barrier  is  due 
to  magnetic  anisotropy,  either  intrinsic  or  shape-induced.  Anisotropy  makes  it 
energetically  possible  for  each  atomic  moment  to  be  aligned  along  one  of  the  “easy” 
axes.  We  consider  a simple  scenario  where  the  total  energy  consists  of  the  leading 
order  term  in  the  anisotropy  energy  and  the  Zeeman  term  (magnetic  field  is  applied 
along  the  z direction) 


E = KiVsin^  0 - M,VH  cos  6, 


(1.26) 


where  V is  the  volume.  The  energy  is  plotted  as  a function  of  9 in  Fig.  1.14. 
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Figure  1.14:  Energy  as  a function  of  the  angle  between  the  2 axis  and  the  magne- 
tization. the  two  minima  are  at  0 = 0 and  at  ^ = tt. 

There  are  two  minima,  one  at  ^ = 0 and  one  at  0 = tt  with  energies  Ei  = 
—MsVH  and  E2  = MgVH  respectively  with  an  energy  barrier  between  them.  At 
thermal  equilibrium,  the  magnetization  will  be  in  the  vicinity  of  these  minima. 
The  peak  value  of  the  energy  is  given  by 


The  nucleation  field  Hn  can  be  obtained  by  minimizing  the  total  energy  (Eq. 


Due  to  thermal  fluctuations  particles  can  jump  from  one  minimum  to  the  other. 
The  number  of  particles  jumping  over  the  barrier  from  minimum  1 to  minimum  2 
per  unit  time  can  be  written  as 


21 


(1.27) 


(1.26)) 


(1.28) 


^12  — Ci2C 


-0(Em-Ex) 


(1.29) 
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where  C12  is  a constant  and  /3  is  related  to  the  temperature  (T)  via  /3  = 1/(A:bT) 
(fcg  is  the  Boltzmann  constant).  Similarly,  the  number  of  particles  jumping  over 
the  barrier  from  minimum  2 to  1 per  unit  time  is 

1^21  (1.30) 

where  C21  is  another  constant.  For  the  case  of  zero  magnetic  field  {H  = 0)  the 
barrier  is  the  same  in  either  directions,  and  these  two  constants  C12,  C21  are  the 
same.  The  relaxation  time  r which  is  the  average  time  the  system  takes  (for 
i/  = 0)  to  jump  from  one  minimum  to  another  is  given  by 

i = (1.31) 

where  /o  is  a constant  of  the  order  of  10^°  s“^  which  can  be  obtained  by  solving 
Brown’s  equation.  The  essence  of  this  simple  calculation  for  the  nucleation  field 
and  the  relaxation  time  is  that  both  can  be  changed  by  changing  the  temperature 
or  the  anisotropy  (and  also  the  magnetic  field  in  the  general  case). 

The  dependence  of  the  energy  of  a ferromagnetic  metal  on  the  direction  of 
magnetic  moment  is  known  as  magnetocrystalline  anisotropy  (MCA)  and  is  one 
of  the  most  important  properties  of  ferromagnetic  materials.  Specially  in  thin 
films  a large  anisotropy  is  seen  to  be  present  at  the  surface.  This  is  now  com- 
monly known  as  surface  anisotropy  and  has  attracted  attention  from  both  theo- 
rists and  experimentalists  [140,  142,  146,  147].  The  total  energy  due  to  surface 
anisotropy  (typically  0. 1-1.0  meV/atom)  is  also  few  orders  of  magnitude  larger 
than  the  total  energy  due  to  bulk  anisotropy  (few  /LxeV/atom)  per  unit  length.  The 
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reduced  symmetry,  the  lower  coordination  number,  and  the  role  of  highly  localized 
surface  states  offer  new  phenomena  which  may  have  applications  in  future  elec- 
tronic devices.  As  an  example  it  has  been  seen  in  photoemission  experiments  [148] 
that  in  a ferromagnetic  thin  film  the  magnetic  moment  is  enhanced  at  the  sur- 
face [148,  149].  Using  SQUID  magnetometry  enhanced  magnetic  moments  for 
Fe  on  Ag(lOO)  with  interface  moments  enhanced  by  as  much  as  29%  have  been 
observed  recently  [150].  The  reduced  coordination  number  at  the  surface  leads 
to  narrowing  of  the  d bands  which  thus  increases  the  density  of  states  at  the 
Fermi  level.  This  is  reflected  in  the  enhancement  of  the  magnetic  moment  at  the 
surface.  Self-consistent  electronic  structure  calculations  also  predict  this  unusual 
feature  [151,  152,  153,  154,  155,  156,  157,  158]. 

Attempts  to  understand  the  origin  of  the  magnetocrystalline  anisotropy  started 
many  years  ago  [159].  More  than  60  years  ago  Van  Vleck  [137]  proposed  that 
the  origin  of  anisotropy  is  a relativistic  effect  due  to  the  coupling  of  the  spin 
of  the  electrons  with  the  angular  momentum.  The  magnetic  dipole  and  other 
possible  anisotropic  exchange  interactions  also  play  roles  in  the  magnetocrys- 
talline anisotropy.  In  transition  metals  the  spin-orbit  coupling  is  smaller  than 
the  crystalline  field  which  makes  the  theoretical  calculations  much  more  difficult. 
Many  theoretical  attempts  using  ab  initio  calculations  [160,  161,  162],  perturba- 
tive techniques  based  on  tight-binding  model  [163,  164,  165],  non-perturbative 
techniques  [166,  167],  etc.  have  been  employed  to  calculate  the  anisotropy  energy 
which  is  several  orders  of  magnitude  smaller  than  other  contributions  to  the  total 
energy  of  a crystal.  All  these  calculations  indicate  that  the  anisotropy  energy  is 
extremely  sensitive  to  the  details  of  the  band  structure  at  the  Fermi  level  - thus 
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correct  estimation  of  the  anisotropy  energy  is  still  a numerical  challenge.  The 
dependence  of  the  MCA  on  the  number  on  layers  also  have  been  studied  both  ex- 
perimentally [140]  and  theoretically  [164,  166,  168,  169,  170].  It  has  been  seen  that 
for  thin  Fe  films  deposited  on  Cu(OOl)  a transition  from  magnetization  perpendic- 
ular to  film  to  magnetization  in-plane  occurs  at  5 ML  [140,  141].  On  the  other 
hand,  no  such  transition  occurs  for  Ni  and  Co  slabs.  At  this  point  it  is  worthwhile 
to  mention  that  in  the  anisotropy  energy  there  are  two  competing  effects  - one 
due  to  spin-orbit  coupling  which  prefers  the  magnetization  to  be  perpendicular 
to  the  plane  and  the  shape  anisotropy  which  prefers  the  magnetization  to  be  in- 
plane. Shape  anisotropy  originates  from  magnetic  dipole-dipole  interactions  and 
for  a monolayer  the  energy  due  to  shape  anisotropy  is  given  by 

■F'shape  ~ 2^^*  ' 

where  the  field  H = N • M*.  N is  known  as  the  demagnetization  factor  which 
depends  on  the  geometry  of  the  sample.  As  an  example  for  thin  films  = Ny  = 
0,  = 1 and  the  shape  anisotropy  energy  is  simply 

= ^^^"cos2  0.  (1.33) 

For  a multilayer  the  generalization  of  Eq.  (1.32)  can  be  found  in  Ref.  [160].  Using 
either  ferromagnetic  resonance  (FMR)  techniques  or  by  torque  measurement,  shape 
anisotropy  can  be  measured  experimentally  [138,  139,  140,  141]. 

For  Fe  films  less  than  five  monolayers  the  anisotropy  due  to  spin-orbit  coupling 
wins  - thus  the  magnetization  prefers  to  be  perpendicular  to  the  plane.  For  more 


33 


than  five  monolayers  the  shape  anisotropy  wins  - thus  the  magnetization  prefers 
to  be  in-plane.  Experimentally  the  measurement  of  the  thickness  of  the  thin  film 
where  this  reorientation  transition  takes  place  depends  on  the  preparation  of  the 
sample.  Theoretically  better  numerical  techniques  are  required  to  explain  the 
experimental  results. 

Landau  and  Lifshitz  [171]  suggested  the  possibility  that  a crystal  placed  in 
an  applied  electric  field  may  experience  a magnetization  which  is  proportional  to 
the  electric  field.  This  magnetoelectric  effect  was  first  observed  experimentally  by 
Astrov  in  a monocrystal  of  antiferromagnetic  Cr203  [172].  It  was  also  observed 
by  Rado  et  al.  that  in  spin-ordered  materials  such  as  Fe304  the  induced  electric 
polarization  is  proportional  to  the  applied  magnetic  field  [173,  174].  The  latter 
effect  is  known  as  magnetically  induced  magnetoelectric  effect  [(ME)^]  and  the 
former  is  known  as  electrically  induced  magnetoelectric  effect  [(ME)^].  To  linear 
order  what  this  means  is  that  the  magnetization  of  the  crystal  is  given  by 

(1.34) 

and  the  polarization  by 

(1.35) 

where  and  H®’'*  are  the  applied  external  electric  and  magnetic  fields  and  oi,P,x 
are  the  proportionality  constants  (summation  over  j is  implied).  Fig.  1.15  shows 
the  result  of  the  experiment  on  magnetite  crystals  by  Rado  et  al.  [174] 

Simple  phenomenological  theory  [173,  174]  as  well  as  microscopic  theory  [175] 
based  on  perturbation  have  been  proposed  to  explain  the  magnetoelectric  phe- 
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Figure  1.15:  Experimental  dependence  of  the  polarization  P on  applied  magnetic 
field  for  different  values  of  the  angle  (ip)  between  the  applied  magnetic  field  and 
the  easy  axis  [174], 
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nomenon  in  spin-ordered  materials.  The  microscopic  origin  of  this  magnetoelectric 
effect  lies  in  the  spin-orbit  coupling  term.  We  can  think  of  the  spatial  part  of  the 
wave-function  as  being  distorted  by  the  applied  electric  field.  Due  to  the  spin-orbit 
coupling  the  spins  are  affected  by  this  distortion  and  hence  they  give  rise  to  an 
induced  magnetization.  In  other  words,  the  electric  field  affects  the  atoms  near  the 
surface  - thus  the  surface  anisotropy  changes  which  again  changes  the  total  energy 
and  the  magnetization  of  the  system. 

Current-driven  excitation  of  magnetic  multilayers  also  have  been  studied  in 
recent  years  [176,  177,  178].  It  has  been  shown  theoretically  that  magnetization  of 
a multilayer  can  be  changed  by  applying  a spin-polarized  current.  Is  it  possible  to 
reverse  the  magnetization  of  a thin  magnetic  film  by  applying  a reasonable  electric 
field?  How  does  the  surface  magnetization  change  with  the  electric  field?  Chapter 
4 addresses  the  above  questions. 

Finally,  Chapter  5 will  briefly  review  the  conclusions  of  the  previous  chapters 
and  provide  suggestions  for  future  work. 


CHAPTER  2 

GIANT  MAGNETORESISTANCE  (GMR)  IN  LATERALLY  CONFINED 

MAGNETIC  MULTILAYERS 

The  importance  of  studying  the  transport  properties  in  magnetic  multilayers 
was  discussed  in  Section  1.1  of  Chapter  1.  The  discussion  was  primarily  focussed 
on  the  new  phenomenon  known  as  the  giant  magnetoresistance  (GMR).  Various 
experimental  results  and  different  theoretical  models  were  also  briefly  outlined.  In 
this  chapter  we  focus  on  GMR  in  laterally  confined  magnetic  multilayers.  In  spite  of 
different  existing  models  - classical  or  quantum  [17,  179]  the  semiclassical  model  of 
Fuchs  and  Sondheimer  [21,  22]  is  quite  successful  in  predicting  the  qualitative  fea- 
tures of  GMR.  Moreover,  one  of  the  strength  of  using  the  Boltzmann  equation  is  to 
get  simple  explanations  of  more  complex  results.  Thus  we  will  use  the  semiclassical 
theory  for  our  calculations  on  laterally  confined  multilayers  in  the  current-in-plane 
(CIP)  geometry.  The  organization  of  this  chapter  is  pedagogical.  We  start  with 
the  semi-classical  Fuchs-Sondheimer’s  derivation  of  conductivity  for  a thin  film. 
We  then  generalize  the  thin-film  result  to  magnetic  multilayers.  Dependencies  of 
GMR  on  different  parameters  are  studied.  The  derivation  for  a single  wire  follows 
next.  To  study  the  effect  of  lateral  confinement  we  compare  the  conductivity  with 
that  of  a film.  Finally  we  generalize  the  single-wire  to  laterally  confined  magnetic 
multilayers.  The  effect  of  lateral  confinement  on  GMR  is  studied  in  detail.  All  our 
results  are  then  summarized  in  the  conclusion. 
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2.1  Semi- Classical  Fuchs-Sondheimer’s  Theory  for  Thin  Films 

Fuchs-Sondheimer  theory  [22]  for  thin  films  is  based  on  the  classical  steady- 
state  Boltzmann  transport  equation: 

V'V./(v,r)-^E-V./(v,r)=  , (2,1) 

' / scatt 

where  /(v,  r)  is  the  distribution  function  for  electrons  of  mass  m at  position  r with 
velocity  v in  presence  of  the  electric  field  E.  Because  of  the  scattering  term,  Eq. 
(2.1)  is  complicated.  For  simplicity  we  consider  the  caise  where  the  scattering  term 
is 

V Acatt  - ’ 

with  (/(v,r))  being  the  spherical  average  of  the  distribution  function,  and  r is  the 
relaxation  time.  This  is  the  simplest  scattering  term  to  represent  elastic  scattering. 
To  solve  Eq.  (2.1)  with  the  right  hand  side  of  Eq.  (2.2),  we  define  ^(v,r)  to  be 
the  deviation  of  the  distribution  function  from  its  equilibrium  value: 

/(v,r)  = /eq(|v|)  + ^(v,r).  (2.3) 

We  next  make  the  ansatz  that  within  linear  response  the  spherical  average  of  the 
distribution  function  is  equal  to  the  equilibrium  distribution  function,  (/(v,r))  = 
/eq(|v|),  i.e.,  the  spherical  average  of  ^(v,  r)  is  zero.  This  ansatz  will  be  tested  later. 
This  approximation  leads  to  a form  known  as  the  relaxation  time  approximation. 
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The  linearized  Boltzmann  equation  then  reduces  to 

V ■ V,(,(v.  r)  - -E  ■ Vv/.<,(|v|)  = - (2.4) 

We  now  consider  the  case  of  a metallic  thin  film  of  thickness  d as  shown  in  Fig. 
2.1.  The  2-axis  is  perpendicular  to  the  plane  of  the  film  and  the  current  is  along 
the  X direction.  In  this  geometry  Eq.  (2.4)  simplifies  to: 


Figure  2.1:  Thin  film  geometry.  Current  is  along  the  x direction. 


dz 


+ eEvx 


g(v,^) 

r 


The  general  solution  of  Eq.  (2.5)  is 


(2.5) 


—erVxE 


-I-  F(v)  exp  ( — , 
V TVJ 


(2.6) 


where  the  energy,  e„  = mv'^/2.  F{v)  is  an  arbitrary  function  of  v to  be  determined 
from  the  boundary  conditions.  At  each  surface  (z  = 0 and  z = d)  the  electrons 
scatter  either  specularly  or  diffusively.  We  assume  that  a fraction  S of  the  electrons 
is  scattered  elastically  at  the  surfaces  and  the  rest  (1-5)  is  scattered  diffusively. 
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5 = 0 corresponds  to  complete  diffusive  scattering  whereas  5 = 1 is  for  perfectly 
specular  scattering.  In  general,  the  specularity  factor  5 can  depend  on  the  direction 
of  the  electrons  but  here  for  simplicity  we  take  it  to  be  a constant.  Depending  on 
the  directions  there  are  two  distribution  functions  g~^  for  electrons  moving  away 
from  the  surface  (with  > 0)  and  g~  for  electrons  coming  into  the  surface  (with 
V2  < 0).  The  boundary  conditions  at  the  z = 0 surface  are 


where  Ci,  C2  are  due  to  diffusive  scattering  at  the  surfaces  and  they  are  indepen- 
dent of  the  direction  of  the  velocities  of  the  electrons.  For  convenience  we  define, 


Combining  Eq.  (2.6)  with  Eq.  (2.9)  the  boundary  conditions  in  Eq.  (2.7)  and  Eq. 
(2.8)  can  be  expressed  as 


From  Eqs.  (2.10)  and  (2.11)  F+  and  F~  can  be  written  in  terms  of  Ci  and  C2. 
Again  at  each  surface  current-in  is  equal  to  current-out  in  the  z-direction.  This 
gives  two  more  equations  to  solve  for  Ci  and  C2-  The  average  of  g^  is  zero  because 
it  is  odd  under  Vx  —Vx-  The  equations  for  C\  and  C2  are  homogeneous  and  the 


^+(v,z  = 0-t-)  = Sg  (v,z  = 0-) -t-Ci, 
g (v,  z = d— ) - Sg^  (v,  z = d-\-)  H-  C2, 


(2.7) 

(2.8) 


(2.9) 
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only  solution  is  Ci  — €2  = 0.  The  distribution  functions  then  are: 
y+(v,2)  = 

1 - 5exp(-d/r|u^|)  ( rlv^l)  }’ 

l-5exp(d/r|u^|)^^^(  r|u^|)}‘ 

By  explicitly  doing  the  spherical  average  of  the  above  distribution  functions  we 
find  that  < g'^(v,  z)  >=  0 and  < ^“(v,  2)  >=  0.  Thus  our  ansatz  is  verified. 

The  current  density  along  the  direction  of  the  electric  field  is  given  by 

Jx{z)  = -2e  j d^vvxg{v,  z).  (2.14) 

The  factor  of  two  accounts  for  the  spin-degeneracy.  The  conductivity  a per  unit 
area  is  obtained  by  averaging  the  current  density 

^ = ^1  dzMz)-  (2.15) 


Restoring  to  polar  coordinates: 

V2  = V cos  9,  Vx  = v sin  6 cos  (f>,  Vy  = v sin  9 sin  (f>, 
d^v  = v^dv  sin  9d9dcj), 
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the  ratio  of  the  conductivity  of  the  film  to  that  of  the  bulk  metal  can  be  derived. 
The  result  is 


^film 

^bulk 


= ( 


cos  6)  cos  6 sin^  6 


1 — exp(— d/L|  cos0|) 

[1  — 5 exp(— d/L|  cos^l)] 


(2.16) 


where  L is  the  mean  free  path  of  the  electrons  and  the  bulk  conductivity  is  defined 
to  be 


ne^L\ 
mvp  j 


(2.17) 


For  perfectly  diffusive  scattering  at  the  surfaces  (S=0),  Eq.  (2.16)  takes  an  ana- 
lytical form  in  terms  of  the  exponential  integral  function  Ei{d/L)  [180] 


^bulk 


(2.18) 


Limiting  values  for  thick  films  and  very  thin  films  can  now  be  obtained  from  Eq. 
(2.16).  For  large  d/L  (thick  films) 


O'film  _ ^ _ 3L(1  — 5) 

<7bulk 


(2.19) 


and  for  thin  films  [d/L  <<  1)  we  get 


C^bulk 


g(l+2S)  (log  ^ + 0.423). 


(2.20) 
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In  the  limit  d/L  » 1 [Eq.  (2.19)]  there  are  two  contributions  to  the  scattering 
rate  - one  due  to  the  background  scattering  and  the  other  due  to  surface  scattering. 
On  the  other  hand,  in  the  limit  d/L  « 1 [Eq.  (2.20)]  the  conductivity  diverges 
(the  resistivity  goes  to  zero)  for  fixed  d when  the  background  scattering  vanishes 
(i.e.  L — >•  oo).  This  is  puzzling  as  one  would  expect  a finite  residual  resistivity 
due  to  the  surface  scattering.  In  quantum  mechanical  models,  this  divergence  of 
the  conductivity  for  large  mean  free  path  disappears  and  the  resistivity  remains 
finite  for  finite  thickness.  This  has  been  explained  in  Section  1.1  of  Chapter  1. 
In  Fig.  2.2  we  plot  the  ratio  of  the  conductivity  of  a film  to  that  of  the  bulk 


Thickness  (in  A°) 


Figure  2.2:  The  ratio  of  conductivity  of  a thin  film  (cTfiim)  to  that  of  the  bulk 
conductivity  ((Xbuik)  as  a function  of  the  thickness  of  the  film.  The  mean  free  path 
of  the  electrons  is  taken  to  be  50  A.  For  specular  scattering  there  is  no  change 
in  the  conductivity  where  as  for  diffusive  scattering  conductivity  increases  with 
thickness  and  finally  saturates  to  the  bulk  value. 
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conductivity  for  three  different  values  of  the  specularity  factors.  For  specular 
scattering  (5  = 1)  there  is  no  change  in  the  conductivity  as  a function  of  thickness 
and  the  conductivity  is  equal  to  the  bulk  value.  But  for  diffusive  scattering  at 
the  surfaces  conductivity  increases  with  thickness  and  finally  saturates  to  the  bulk 
value. 


Calculation  of  the  conductivity  and  the  giant  magnetoresistance  (GMR)  of  a 
magnetic  multilayer  structure  is  due  to  Camley  and  Barnas  [23,  24].  Their  theory 
is  the  generalization  of  Fuchs-Sondheimer  theory  for  thin  films.  The  Boltzmann 
equation  [Eq.  (2.1)]  is  same  as  before  with  the  only  difference  that  the  distribution 
functions  and  the  parameters  carry  an  extra  index  n for  the  layers.  A N layer 
system  which  is  alternating  layers  of  ferromagnets  (FM)  and  paramagnets  (PM)  is 
shown  schematically  in  Fig.  2.3.  The  general  solution  for  the  distribution  function 
is  thus 


The  functional  form  of  E„(v)  is  determined  from  the  boundary  conditions.  We 
proceed  in  the  same  fashion  as  we  did  for  the  film.  Let  g^{\,z)  and  g~(\,z)  be 
the  distribution  functions  of  the  electrons  with  > 0 and  with  < 0 respectively. 
At  each  interface  the  electrons  undergo  specular  or  diffuse  scattering.  For  specular 
scattering,  the  reflection  and  transmission  coefficients  enter  into  the  boundary 


2.2  Generalization  to  Magnetic  Multilayers 


(2.21) 
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interface: 

^N=  d]+d2+d3+dN 
Z 3=  dj"*"d2'*’d3 

Z2=  di  +d2 

Zi=di 

Zo=  0 

Figure  2.3:  Schematic  diagram  of  a N-layer  system.  Thickness  of  the  n-th  layer  is 

dfi- 

conditions.  The  boundary  conditions  at  the  nth  interface  are 


9n  (v,  Z-^)  = SnTnQn-l  (v,  ) + S^R^g^  (v,  Z+ ) 

1 + C„, 

(2.22) 

9n-l  (v,  Z~)  = SnTng-  (v,  2+)  + SnRng^-l  (v,  Z~) 

1 + c„. 

(2.23) 

where  and  T„  are  reflection  and  transmission  coefficients  at  the  nth  interface, 
and  Sn  is  the  specularity  factor.  C„  is  a constant  due  to  the  diffuse  scattering  at  the 
interface.  At  each  interface  the  sum  of  the  reflection  and  transmission  coefficients 
is  unity,  i.e.  Rn  T„  = 1.  We  have  assumed  T,  R,  and  C depend  only  on  the 
interface  - not  on  how  the  electrons  enter  or  exit  the  interface.  Also  to  mention 
that  the  transmission  coefficient  is  zero  (thus  R=  1)  for  the  first  (n  = 0)  and  the 
la^t  surface  (n  = AT). 
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We  now  combine  the  boundary  conditions  in  Eqs.  (2.22),  (2.23),  and  the  solu- 
tion of  the  Boltzmann  equation  [Eq.  (2.21)]  within  each  layer.  For  convenience  we 
choose  the  following  convention: 

9n  (v, z)  = gl{v)  + F+  (v)  exp  , (2.24) 

9n  (v,  z)  = gl(v)  + F~  (v)  exp  j , (2.25) 

where  we  have  introduced  the  simplifying  notation: 

i/n(^)  = -eEvxTn  ^ • (^  = |v|)  (2.26) 

This  allows  us  to  express  the  boundary  conditions  solely  in  terms  of  the  F’s  and 
C’s: 


5,!'(v,4)  = 5°(u)  + F+(v),  (2.27) 

9n-i  (v,  Z~)  = ^“(v)  -I-  (v)  exp  ) , (2.28) 

V Pn-l  / 

9n  (v,  2;T)  = ^°(v)  + F~  (v)  exp  , (2.29) 

= ^°(^^)  + (v) . (2.30) 


Substituting  Eqs.  (2.27)  and  (2.30)  in  Eqs.  (2.22)  and  (2.23)  yields: 

SnTngl_M  - {I  - SnRn)g»  - Cr.  = F+ (v) 

- 5„T„F+_i  (v)  exp  ) 
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- S^RnF^  (v)exp 

SnTMv)  - {1  - SMgl_,{v)  - C,,  = F„t,(v) 

- SnRr,F^_i  (v)  exp 


5„T„F„  (v)exp 


There  are  2N  of  these  equations  - one  for  each  F^.  We  go  through  them  sequen- 
tially and  eliminate  the  F^  one  at  a time.  At  each  step,  we  move  terms  to  the 


or  more  of  the  exponential  functions  with  |u^|.  On  the  right-hand-side  are  terms 
which  are  linear  in  the  They  also  may  be  multiplied  by  one  or  more  of  the 
exponential  functions  with  |uz|.  This  means  that  the  final  F^  and  each  of  the 
others  have  the  form 


This  still  leaves  N + 1 number  of  C’s.  The  remaining  N + 1 equations  are  that 
the  current-in  is  equal  to  the  current-out  in  the  z-direction  at  each  interface.  The 
average  of  the  ^°’s  is  zero  because  they  are  odd  under  Vx  -Vx-  The  equations 
for  the  C's  are  homogeneous,  and  baring  any  special  coincidence,  the  only  solution 
is  that  all  the  C’s  are  zero. 

The  2N  equations  we  must  now  solve  are: 


left-hand-side  containing  and  the  C’s.  These  terms  may  be  multiplied  by  one 


F^  = 


E {g°'s)  X (exp's)  -f  E (C's)  x (exp’s) 
E (exp's) 


(2.33) 
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- (v)exp  (--^y  (2.34) 

V Pzl-Tn/ 

- (1  - = F„ti(v)-5„T„F-(v)expf--^) 

V \Vz\rnJ 

- 5„i?„F+_i  (v)exp  (-  ^ . (2.35) 

We  now  simplify  the  above  equations.  Each  of  the  g°'s  on  the  left-hand-side  has 
the  common  factor  —eEvx{—dfeq/dev)  which  we  can  divide  out.  Defining 


F(v)  = -eB«. 


(2.36) 


Eqs.  (2.34)  and  (2.35)  now  read 


- (1  - 5„F„)/„  = F+(v)-5„r„F„ti(v)expf-^%:^ 


- 5„F„F„  (v)exp  ^ 
V Vz  T„ 


(2.37) 


5„T„/„  - (1  - 5„F„)f„_x  = F+_x(v)-5„T„F-(v)expf-^) 

V l^zPn/ 


- 5„F„F+_x  (v)  exp 


dn-i  ' 

kzk„-i. 


(2.38) 


where  /„  = vpTn  is  the  mean  free  path  in  the  nth  layer.  The  angular  dependence 
in  these  equations  only  comes  in  through  \vz\.  Letting  u=  | cos^|,  we  get 


SnTrdn-l  ~ (1  ~ SnRn)ln 


SnTJn  - (1  - 5„F„)f„_ 


Ki'^)  - -S„T„F+_x(u)exp 

V L-iu/ 

S„RnF~{u)exp  > (2-39) 

Fn-i{u)  - 5„r„F-(u)exp 


- 5„i?„F+_i(u)exp 


^n— 1 

In-lU 
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(2.40) 


The  above  equations  have  to  solved  for  u between  zero  and  one. 

Averaging  the  current  over  the  z-direction  the  conductivity  can  be  expressed 
as 

a = ^ f dzj^{z)  = N{EF)e^VF  Y, 

+ JV(E^)eV|; /‘f  i (l  - u^)  U {l  - exp  (-Aj  | 

X (u)  + K i^)}  ■ (2.41) 

where  N{Ef)  is  the  density  of  states  at  the  Fermi  surface.  For  this  multi-layer 
geometry  we  compute  the  conductivity  in  two  cases:  when  the  ferromagnet’s  mag- 
netizations are  parallel,  ct^t,  and  when  they  are  antiparallel,  ctaf-  The  magnetore- 
sistance is  defined  as  the  ratio  of  the  change  in  the  conductivity  from  parallel  to 
antiparallel  alignments  divided  by  the  parallel  conductivity: 


GMR  = 1 - — . (2.42) 

(Tp 

In  Fig.  2.4  we  plot  the  giant  magnetoresistance  of  a three-layer  structure  as 
a function  of  thickness  of  the  spacer  layer  for  three  cases:  (a)  GMR  due  to  bulk 
scattering,  (b)  GMR  due  to  surface  scattering,  and  (c)  GMR  due  to  both  bulk  and 
surface  scattering.  For  bulk  scattering  [case  (a)]  the  GMR  decreases  monotonically 
whereas  for  surface  scattering  [case  (b)]  the  GMR  initially  increases  and  then 
decrecises  with  the  thickness  of  the  spacer  layer.  This  can  be  understood  from  the 
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Figure  2.4:  GMR  of  a three-layer  structure  (Co/Cu/Co)  as  a function  of  the  thick- 
ness of  the  Cu  layer  (tcu)  for  (a)  bulk  scattering,  (b)  surface  scattering,  and  (c) 
both  bulk  and  surface  scattering.  The  thickness  of  the  Co  layers  are  kept  fixed  at 
8 A.  For  cases  (a)  and  (c)  the  mean  free  paths  are  chosen  to  be  that  of  Co/Cu/Co, 
— 55  A,  = 10  A,  and  = 226  A [181].  For  case  (b)  in  which  the  surface 
scattering  contributes  to  the  GMR  we  use  L'’"  = L-*-  = 17  A and  = 226  A. 
For  cases  (b)  and  (c)  the  transmission  coefficients  are  taken  to  be  = 0.8,  = 

0.4  [182].  In  case  of  (a)  which  is  for  the  bulk  scattering  we  take  T'’"  = T-*-  = 0.6. 
We  choose  the  interfaces  to  have  the  same  specularity  factor  S = 0.9. 
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classical  model  we  derived  in  Chapter  1.  Following  Eq.  (1.10)  we  see  that  for  bulk 
scattering  in  the  limit  of  zero  thickness  of  the  spacer  layer  the  GMR  is  finite  and  is 
equal  to  the  GMR  obtained  from  Eq.  (1.5).  The  decrease  in  GMR  with  increase  in 
thickness  of  the  spacer  layer  also  follows  from  Eq.  (1.10).  The  maximum  value  of 
the  GMR  for  case  (a)  primarily  depends  on  the  mean  free  paths  of  the  ferromagnet. 
The  general  case  of  both  bulk  and  surface  scattering  [case  (c)j  is  a combination  of 
case  (a)  and  case  (b).  The  decrease  of  GMR  with  the  thickness  of  the  spacer  layer 
is  in  agreement  qualitatively  with  the  experimental  results  [9,  10,  11]. 

GMR  as  a function  of  thickness  of  the  ferromagnet  for  the  same  cases  as  in 
Fig.  2.4  is  plotted  in  Fig.  2.5.  In  all  cases  GMR  initially  increases  and  then 
decreases  with  increase  in  the  thickness  of  the  ferromagnetic  layers.  Also  note  that 
the  contribution  to  GMR  due  to  bulk  scattering  is  more  than  the  contribution 
from  the  interface  scattering. 

The  eflFect  of  the  surface  roughness  (5)  on  the  GMR  is  shown  in  Fig.  2.6.  With 
entirely  diffusive  (S  = 0)  surfaces  and  interfaces  the  multilayers  are  independent 
of  each  other  which  means  that  there  is  no  GMR.  With  a decrease  in  the  interface 
roughness,  the  GMR  increases  and  finally  reaches  a maximum  for  perfectly  specular 
surfaces. 


2.3  Thin  wires  - Theory  and  Results 

In  the  case  of  thin  wires  the  problem  is  two-dimensional  and  the  Boltzmann 
equation  takes  the  form  of  a partial  differential  equation  in  the  space  variables  [183]. 
With  completely  diffusive  surfaces  results  for  wires  with  square  cross-section  were 
obtained  by  MacDonald  and  Sarginson  [184]  and  for  wires  with  circular  cross- 
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Figure  2.5:  GMR  of  a three-layer  structure  (Co/Cu/Co)  as  a function  of  the  thick- 
ness of  the  Co  layer  (tco)  for  (a)  bulk  scattering,  (b)  surface  scattering,  and  (c) 
both  bulk  and  surface  scattering.  The  thickness  of  the  Cu  layer  is  kept  fixed  at 
10  A.  In  all  cases  GMR  initially  increa.ses  and  then  decreases  with  increase  in  the 
thickness  of  the  ferromagnet.  Note  that  GMR  due  to  bulk  scattering  is  more  than 
the  surface  scattering.  The  parameters  chosen  are  the  same  as  in  Fig.  2.4. 
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Figure  2.6:  GMR  is  plotted  as  a function  of  the  specularity  factor  for  different 
thicknesses  (t)  of  the  layers.  GMR  increases  with  decrease  in  the  interface  rough- 
ness. For  convenience  the  thickness  of  the  Co  and  Cu  layers  is  taken  to  be  the 
same.  We  use  the  same  parameters  as  in  Fig.  2.4. 
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section  by  Dingle  [185].  In  this  section  we  first  present  analytical  results  for  the 
conductivity  of  a thin  wire  with  rectangular  cross-section  but  with  diflPusive  sur- 
faces. We  then  generalize  to  wires  with  specular  and  diffusive  surfaces.  In  the 
latter  case  we  solve  for  the  conductivity  numerically  using  an  iterative  procedure. 

We  consider  a wire  of  width  w and  thickness  t.  For  convenience,  the  electric 
field  is  taken  to  be  along  the  2 direction.  Similar  to  Eq.  (2.5)  the  Boltzmann 
equation  in  this  geometry  under  the  relajcation  time  approximation  reads 


_ dg{v,x,y)  , dg{v,x,y)  , g{v,x,y)  ^ ( 5/eq' 

Vx h 1 = —ebjV,  ' 


OX 


dy 


(2.43) 


The  general  solution  of  Eq.  (2.43)  can  be  written  as 


g{v,  X,  y)  = -eEv^  ^ ^ f — ~ 

V / L V TVxJ  \Vy 


y ^ 

Vx, 


(2.44) 


or  equally  well, 


^(v,  X,  y)  = -eEv^ 


df, 


eq 


de„ 


1 — exp  I — 


TV,j 


\Vx  Vy^ 


(2.45) 


The  functions  ^ and  $ are  arbitrary  and  have  to  be  determined  from  the  boundary 
conditions.  We  assume  scattering  is  completely  diffusive  at  the  boundaries  - so  the 
distribution  function  g of  the  electrons  leaving  the  boundaries  must  be  independent 
of  the  ratio  Vx  : Vy  : at  x = 0,  w,  and  y = 0,t.  Specifically  g is  equal  to  zero 

under  these  conditions.  The  complete  distribution  function  can  be  represented  by 
the  following  set  of  equations  in  velocity  space: 
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First  quadrant:  Vx  > 0,  Uy  > 0, 


g > 0,  Uy  > 0,  x,  y < = -eEv^ 


1 — exp 


TV. 


y / \ 


g > 0,  Uy  > 0,  X,  y > ^x^  = -eEv^ 

1) 

TVxJ 


1 — exp 


Second  quadrant:  u*  > 0,  Uy  < 0, 


y Uj;>o,Uy<o,x< — -{t-y),y\  = -e£;u^(- 


de„ 


y Uj  > 0,Uy  < 0,  X > {t-y),y]  = -eEv, 


1-exp  r — , 

\ TVx ) \ 

^ dfeq\ 


X 


1 — exp 


d€v  j 
't  - 2/' 


TV,, 


Third  quadrant:  Uj  < 0,  Uy  < 0, 


y ( Uj,  < 0,Uy  < 0,  (it;  - x)  < — (t  - y),y  ) = -eEv 


df. 


eq 


X J 

^ /w  — X 

1 — exp  

\ TVx 


g [vx  < 0,Vy  < 0,{w  - x)  > —{t  - y),y]  = -eEv^ 


1 — exp 


df. 


eq 


dey  J 
't  - y'' 


TV,, 


(2.46) 

(2.47) 


(2.48) 


(2.49) 


, (2.50) 


X 


. (2.61) 


55 


Fourth  quadrant:  Vx  < 0,Vy  > 0, 


g (vx  < 0,Vy  > 0,x,y  < -^{w  - x)^  = -eEv^ 


X 


eq 


d€y 


1 — exp 


TV,, 


g (vx  < 0,Vy  > 0,x,y  > -^{w  - x)^  — -eEv^ 


df. 


eq 


de„ 


1 — exp 


— x\ 

\ TVx  ) 


(2.52) 


(2.53) 


Due  to  symmetry  each  quadrant  contributes  equally  to  the  total  current  J^.  Using 
Eq.  (2.15)  and  Eqs.  (2.46  - 2.53)  the  conductivity  is  calculated  and  the  result  is, 


/ d^sin^^cos^^  / 

.TTW  I Jo  Jo 


+ 


+ 


ftein  ^(w/t) 


/6L^\  W2  /-t— 

[id  I Jo 

f— ) r^^de  sin^ecos^e  _ e-^/Lsm0 

\TTWt J Jo  L 


f6L^\  W2 
V TTWt  ) Jo 


ddsin^dcos^ 


in0  J 


(2.54) 


Different  limiting  values  of  Eq.  (2.54)  can  be  obtained.  In  the  limit  {L  « w,t) 
we  obtain 

^ « 1 _ (3 WL  ^ L\  /£W^\ 

\8/  \w  t)  \57r/  \ wtj 


(2.55) 
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In  the  other  limit  i.e.  L » w,t  we  get 


•^wire  ^ ^ 3m  ^ ^ \/w^  + ”1"  + m \ 

(Tfilm  \8L/  yy/u)^-\-t2_fj  V8L/  \^/w^  + f^  — wj 


(2.56) 


where  7 is  the  euler  constant. 

For  the  case  of  a wire  with  square  cross-section  (m  = t),  Eqs.  (2.55)  and  (2.56) 
can  be  simplified. 


^wire 


+ {L«t) 


(2.57) 


and  in  the  other  limit  {L  » t) 


^wire 


t 

1.115-. 

L/ 


(2.58) 


Comparing  Eq.  (2.57)  with  the  results  of  a thin  film  with  5 = 0 in  Eq.  (2.19) 
and  Eq.  (2.58)  with  Eq.  (2.20)  the  much  more  rapid  decrease  of  the  conductivity 
in  the  case  of  a wire  is  evident.  The  result  for  the  conductivity  of  a wire  with 
rectangular  cross-section  but  diffusive  surfaces  is  shown  in  Fig.  2.7.  The  ratio  of 
conductivities  is  compared  with  that  of  a thin  film. 
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Figure  2.7:  Comparison  of  a thin  square  wire  result  with  that  of  a thin  film.  The 
ratio  of  the  conductivities  (cr/abuik)  is  plotted  as  a function  of  the  ratio  of  the 
mean  free  path  (L)  and  the  width  (a)  (thickness  in  case  of  the  film).  Due  to 
lateral  confinement  cr/cTbvOk  decreases  faster  in  case  of  wire  than  that  of  the  film. 
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We  numerically  solve  the  general  case,  i.e.  when  scattering  at  the  surfaces  is 
both  specular  and  diffusive.  We  go  back  to  the  linearized  Boltzmann  equation  in 
Eq.  (2.1).  When  the  electric  field  is  zero,  Eq.  (2.4)  becomes  a homogeneous 
equation  with  the  expected  solution  being  g{'v,  r)  = 0.  For  non-zero  electric 
field,  g(v,r)  is  proportional  to  (eE/m)Vv/eq(|v|)  = eE  • va/eq(|v|)/ae„.  Thus, 
g{-v,r)  = -g{v,T),  irrespective  of  the  boundary  conditions,  and  the  ansatz  that 
the  spherical  average  of  g{v,  r)  is  zero  is  justified.  We  write  the  general  solution 
in  a slightly  different  form  than  before  [186,  187] 


^(v,r)  = ^(v,rB)e  I""  + erE  • v (1  - (2.59) 


where  is  a point  on  the  boundary.  Eq.  (2.59)  implies  that  to  find  the  distribution 
function  at  position  r with  velocity  v,  we  proceed  from  r backwards  along  v until 
we  reach  a point  at  the  boundary.  The  distribution  function  at  the  boundary, 
^(v.tb),  is  determined  by  the  boundary  conditions.  We  recover  the  usual  bulk 
value  if  we  go  far  away  from  the  boundary  points.  Also  we  can  see  from  Eq.  (2.59) 
that  the  electrons  lose  their  momentum  as  they  diflFuse  into  the  medium  and  the 
characteristic  length  scale  for  this  is  just  the  mean  free  path,  rvp,  where  vp  is  the 
Fermi  velocity. 

We  now  examine  the  boundary  conditions.  Electrons  are  scattered  either  spec- 
ularly or  diffusively  at  each  surface.  For  a thin  wire  there  are  four  surfaces.  A 
schematic  diagram  of  a rectangular  wire  is  shown  in  Fig.  2.8.  Let  j be  an  index 
for  the  surface.  The  boundary  condition  at  the  j-th  surface  can  be  expressed  as: 


5r(v,rj)  = 5^f(-v,rj), 


(2.60) 
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where  rj  is  a point  on  the  j-th  surface,  and  the  superscripts  out,  in  correspond  to 
electrons  going  out  from  the  boundary  or  coming  in  to  the  boundary.  We  iteratively 
solve  the  linear  set  of  equations  for  the  distribution  functions,  5'°“*(v,rB),  which 
are  non-uniform  along  the  x and  y directions  of  the  surfaces.  Prom  Eq.  (2.4)  we 
observe  that  the  r)’s  with  different  velocities  projected  along  the  z direction 

are  decoupled.  This  allows  us  to  discretize  the  r)  according  to  cos^  = '^z/\v\ 

and  solve  each  separately.  For  each  cosO,  at  the  surfaces  carry  two  more 

indices:  an  angle  (f>  and  a position  r^.  The  p°“*(v,r)’s  for  different  4>  and  Tj  are 
related  via  Eq.  (2.6)  and  the  boundary  conditions  given  in  Eq.  (2.60).  As  an 


z 


-e Width  ^ 


Figure  2.8:  Schematic  diagram  of  a thin  wire.  The  four  surfaces  are  labeled  as  1, 
2,  3,  and  4.  Current  is  in  plane  along  the  z direction.  The  distribution  function 
for  outgoing  electrons  at  point  r*,  on  surface  1 depend  on  all  possible  incoming 
electrons  from  points  on  the  other  surfaces.  As  an  example  the  contribution  from 
point  Fj  on  surface  4 is  shown  schematically.  The  precise  relationship  is  described 
in  the  text  [see  Eq.  (2.61)]. 


example  the  relation  between  the  distribution  functions  of  the  outgoing  electrons 
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at  a surface  1 is: 


^r‘(cos  e,  (f),  Tk)  = Sg^\cos  d,  27t  - cf),  x (2.61) 


where  di  is  the  path  length  projected  onto  the  a;  - y plane  from  one  boundary  to 
another.  Initial  guesses  for  yJ“*(cos6>,  0,ri)  are  taken  from  the  nearby  yJ“‘[cos6>  - 
5(cos0),0,  Fj].  The  calculations  are  converged  to  within  1%. 

Once  we  know  the  distribution  functions  at  the  boundaries,  by  using  Eqs. 


momentum  v at  any  point  r.  The  conductivity  is  then  calculated  by  averaging  the 
current  density  over  an  area  A: 


The  factor  of  two  accounts  for  the  spin-degeneracy. 

In  Fig.  2.9  we  plot  the  ratio  of  the  conductivity  of  the  wire  to  that  of  the  bulk 
as  a function  of  the  ratio  of  the  mean  free  path  to  the  width  of  the  wire.  Similar 
to  thin  films  the  ratio  <T^re/o’buik  decreases  with  the  ratio  L/width,  where  L is  the 
mean  free  path. 

To  show  the  effect  of  lateral  confinement  we  plot  in  Fig.  2.10  the  ratio  of 
conductivities  cr^^^re/crbuik  as  a function  of  the  aspect  ratio  (width/thickness)  for 
different  values  of  the  specularity  factor.  For  perfect  specular  scattering  the  ratio 
is  unity.  In  all  other  cases  the  ratio  increases  with  decrease  in  the  aspect  ratio. 
For  very  large  values  of  the  aspect  ratio  the  conductivity  of  the  wire  saturates  to 
the  bulk  conductivity. 


(2.3)  and  (2.6)  we  can  determine  the  distribution  function  of  the  electrons  with 


(2.62) 
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Figure  2.9:  Conductivity  of  a thin  wire  (uvrire)  with  respect  to  the  bulk  conductivity 
(c’’buik)  is  plotted  as  a function  of  the  ratio  of  the  mean  free  path  (L)  and  the  width 
of  the  wire  for  three  different  values  of  the  specularity  factor  (S).  The  aspect  ratio 
(thickness/width)  is  kept  fixed  to  unity  in  all  three  cases.  Similar  to  thin  films  the 
ratio  of  the  conductivities  decreases  with  increase  in  L/Width. 
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Figure  2.10:  cr^re/c^buik  as  a function  of  the  aspect  ratio  (width/thickness).  The 
ratio  of  the  conductivities  increases  with  increase  in  the  the  aspect  ratio  for  different 
values  of  the  specularity  factor  S {S  ^ 1).  For  5 = 1 the  conductivity  of  the  wire 
is  equal  to  the  bulk  conductivity.  The  thickness  is  kept  fixed  to  10  A in  all  cases 
and  the  mean  free  path  is  chosen  to  be  50  A. 
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2.4  Generalization  to  Laterally  Confined  Magnetic  Multilayers 

The  geometry  of  our  problem  is  shown  in  Fig.  2.11.  We  compute  the  current-in- 
plane conductivity  (CIP)  for  a three-wire  structure  stacked  along  the  ^-direction. 
Both  the  current  and  the  electric  field  are  in  the  z direction.  We  have  two  identical 
ferromagnetic  materials  (FM)  and  a paramagnet  (PM).  For  this  three  wire  geome- 
try we  compute  the  conductivity  in  the  case  when  the  ferromagnet’s  magnetizations 
are  parallel,  ap,  and  when  they  are  antiparallel,  (Taf-  The  magnetoresistance  is 
then  calculated  from  Eq.  (2.42). 

For  the  laterally  confined  magnetic  multilayers  the  solution  of  the  Boltzmann 
equation  carries  an  extra  layer  index  n.  The  generalized  version  of  Eq.  (2.59)  reads 


^n<(v,r)  = gns{v,TB)e  !'■ 


+ 


er„gE  • V 


(2.63) 


where  is  a point  on  the  boundary  or  interface  and  s refers  to  the  spin  of  the 
electrons.  We  now  consider  the  boundary  conditions.  At  each  interface  the  elec- 
trons undergo  either  specular  or  diffuse  scattering.  For  the  nth  interface,  which 
is  between  wires  n and  n -I- 1,  we  define  the  probability  of  spin  s electrons  being 
diffusively  scattered  as  (1  — Sns),  where  Sns  is  the  spin  dependent  specularity  fac- 
tor. The  probabilities  for  being  specularly  refiected  and  transmitted  are  SnsRns 
and  SnsTns,  respectively.  The  sum  of  Rns  and  is  one.  The  angular  depen- 

dence of  the  surface  scattering  parameter,  [188,  189]  the  reflection  coefficients,  and 
the  transmission  coefficients  [30,  182]  have  been  studied  before.  For  simplicity  in 
our  calculations  we  treat  those  as  angle  independent.  With  these  definitions,  the 


64 


Width 


Figure  2.11;  Schematic  diagram  of  the  three  wire  structure.  Wires  1 and  3 are 
the  ferromagnets  of  thickness  and  wire  2 is  a paramagnet  of  thickness  tpM 
stacked  along  the  y direction.  Current  is  in  plane  (CIP)  along  the  z direction.  S 
represents  the  spin  independent  specularity  factors  at  the  sides  of  the  ferromagnet, 
and  the  paramagnet.  Considering  all  possible  incoming  electrons  from  points  on 
the  edges  and  interfaces,  e.g.,  points  r,-  and  r*,  the  distribution  functions  for  the 
outgoing  electrons  at  point  Tj,  g{cos  9,  (j),  tj)  are  determined.  The  exact  relationship 
is  described  in  the  text  [see  Eq.  (2.66)]. 
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boundary  conditions  at  the  nth  interface  can  be  expressed  as: 


where  y„  is  the  position  of  the  interface,  and  the  superscripts  out,  in  correspond 
to  electrons  going  out  from  the  boundary  or  coming  in  to  the  boundary.  Similar 
equations  are  satisfied  at  the  sides  of  the  layers  except  that  there  is  no  transmission 
{T  = ^,R  = 1).  Fig.  2.12  illustrates  one  of  the  boundary  conditions  [Eq.  (2.64)]. 


Figure  2.12:  Pictorial  representation  of  Eq.  (2.64).  Only  the  outgoing  and  incom- 
ing distribution  functions  of  the  electrons  at  the  interface  of  (n  - l)th  and  nth 
layers  are  shown.  The  reflected  part  of  the  incoming  electrons  in  the  nth  layer, 
a;,  y„),  and  the  transmitted  part  of  the  incoming  electrons  in  the  (n  - l)th 
layer,  gn"Lig(v,  x,  yn)  contribute  to  the  outgoing  distribution  function  g^\v,  x,  ?/„). 
The  specularity  factor  of  the  interface  is  5„j. 

We  use  an  iterative  procedure  to  compute  the  distribution  functions,  ^^‘(v,  tb), 
which  are  non-uniform  along  the  x and  y directions  of  the  interfaces  and  edges  (see 
Fig.  2.11).  Our  numerical  procedure  is  similar  to  that  applied  to  calculate  the  con- 


n-1 


g„  (v,x,y„) 


n 
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ductivity  for  a thin  wire  in  Section  2.3.  However  we  briefly  describe  our  numerical 
procedure  below. 

From  Eq.  (2.4),  we  see  that  the  ^^*(v,r)’s  with  different  velocities  projected 
along  the  z direction  are  decoupled.  This  allows  us  to  discretize  the  g^\v,r) 
according  to  cos 9 = Vz/\v\  and  solve  each  separately.  For  each  cos9,  g^*{v,T)  at 
the  edges  and  the  interfaces  carry  two  more  indices:  an  angle  <p  and  a position  rj. 
The  i?^‘(v,r)’s  for  different  (f)  and  rj  are  related  via  Eq.  (2.6)  and  the  boundary 
conditions  given  in  Eqs.  (2.64)  and  (2.65).  To  illustrate  this  we  consider  the 
relations  between  the  distribution  functions  of  the  outgoing  electrons  at  the  first 
interface  in  Fig.  2.11: 

g^\cos  9,  <l>,  rj)  = Si,Ri,g^*(cos  9,  27t  - 

+ (2.66) 

Here  di,  d2  are  the  path  lengths  projected  onto  the  x — y plane  from  one  boundary 
to  another.  Initial  guesses  for  5^‘(cos  9, 4>,  Vi)  are  taken  from  the  nearby  g^J[cos  9- 
5{cos9),  <f),  Fj].  The  calculations  are  converged  to  within  1%  with  the  total  number 
of  divisions  for  cos9,  r*,  and  </»  chosen  as  N^ose  = 100,  = 500,  and  Ni  = 400, 

respectively. 

From  the  distribution  functions  at  the  boundaries  using  Eqs.  (2.3)  and  (2.63) 
we  can  determine  the  distribution  function  of  the  electrons  with  momentum  v 
at  any  point  r.  Following  Eq.  (2.15)  we  then  calculate  the  conductivities  both 
for  parallel  and  antiparallel  alignments  of  the  magnetizations.  Finally,  the  giant 
magnetoresistance  is  obtained  from  Eq.  (2.42). 
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We  now  consider  the  different  input  parameters  of  our  problem.  They  are; 
the  spin  dependent  mean  free  paths  in  the  ferromagnet,  the  mean  free  path 
in  the  paramagnet,  the  spin-dependent  transmission  coefficients,  the 

specularity  factors,  S,  and  the  thicknesses  of  the  layers,  tpM  and  tpu-  We  choose 
the  mean  free  paths  to  be  those  determined  in  an  experiment  on  a Co/Cu/Co 
structure:  [181]  = 55  A,  = 10  A,  and  = 226  A.  The  transmission 

coefficients  are  obtained  from  an  average  of  the  transmission  coefficients  calculated 
by  Stiles:  [182]  = 0.8  and  = 0.4.  We  take  the  specularity  factors  at  all 

interfaces  and  sides  to  be  the  same:  S = 0.9.  Finally,  the  thicknesses  are  chosen 
to  be  those  of  Cu  when  the  Co  layers  are  antiferromagnetically  coupled  at  zero 
field  [9].  For  simplicity  we  take  the  ferromagnets  to  have  the  same  thicknesses  as 
the  paramagnets:  tpM  = ^pm- 

GMR  is  due  to  both  spin  anisotropies  in  the  bulk  and  surface  scattering.  To 
understand  the  role  of  bulk  and  surface  scattering  in  GMR  we  consider  the  limiting 
cases  i.e.,  when  the  GMR  is  due  to  only  bulk  scattering  anisotropies  or  only  surface 
scattering  anisotropies.  To  do  this,  we  consider  two  special  cases:  (a)  when  the 
transmission  coefficients  are  equal:  = 0.6  and  (b)  when  the  bulk  mean 

free  paths  in  the  ferromagnet  are  equal:  = 17  A.  In  the  following  these 

are  referred  to  as  the  (a)  bulk  scattering  case  and  (b)  surface  scattering  case. 

In  Fig.  2.13  we  plot  the  giant  magnetoresistance  as  a function  of  width  for  the 
three  cases:  (a)  GMR  due  to  bulk  scattering  anisotropies,  (b)  GMR  due  to  surface 
scattering  anisotropies,  and  (c)  GMR  due  to  both  bulk  and  surface  scattering 
anisotropies.  In  all  cases  the  giant  magnetoresistance  decreases  with  decrease  in 
the  width.  The  decrease  of  GMR  with  lateral  confinement  can  be  understood  by 
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considering  a film,  i.e.,  a wire  with  infinite  width.  Diflfusive  scattering  at  the  sides 
of  the  wire  reduces  the  conductivity  and  mean  free  paths  of  the  wire  compared 
to  the  film.  This  reduction  in  the  conductivity  in  going  from  a film  to  a wire 
should  be  comparable  to  the  reduction  in  going  from  a bulk  system  to  a film  with 
a thickness  equal  to  w.  Following  Eq.  (2.16)  the  conductivity  of  such  a film  is 


where  the  bulk  conductivity,  cTbxUk  is  defined  in  Eq.  (2.17).  We  obtain  an  effective 
mean  free  path,  L^e,  by  replacing  the  mean  free  path  L in  Eq.  (2.17)  with  L^s 


calculation  with  infinite  width  wires.  The  results  of  such  an  effective  mean  free 
path  multilayer  calculation  (similar  to  Section  2.4)  are  plotted  as  the  solid  lines 
in  Fig.  2.13.  This  approximation  is  in  good  agreement  with  the  exact  results 
(symbols),  and  hence  we  conclude  that  the  reduction  in  the  GMR  is  due  to  a 
decrease  in  the  effective  mean  free  path  within  the  layers. 

Although  the  GMR  is  reduced  in  all  the  cases  shown  in  Fig.  2.13,  the  actual 
GMR  vs.  width  curves  are  different.  In  particular,  there  are  different  length  scales 
at  which  the  GMR  is  reduced.  The  width  at  which  the  GMR  is  half  its  infinite 
width  value,  GMR(wire)/GMR(film)  = 1/2,  is  defined  as  the  half-width.  In  Fig. 
2.14  we  have  rescaled  the  GMR  vs.  width  curves  shown  in  Fig.  2.13  by  GMR(film) 
and  the  half-width.  All  the  points  fall  close  to  a single  curve.  This  means  that  the 
half-width  and  the  GMR  for  infinite  width  wires,  GMR(film),  determine  the  GMR 


(2.67) 


such  that  (Tbuik  = o’fUm-  The  effective  mean  free  paths  can  be  used  in  a multilayer 


69 


Figure  2.13:  Giant  magnetoresistance  of  a three  wire  structure  as  a function  of 
width  for  (a)  bulk  scattering,  (b)  surface  scattering  and  (c)  both  bulk  and  surface 
scattering.  The  symbols  refer  to  different  film  thicknesses,  tpM  = tpM  = 8 A (circle), 
20  A (box),  and  30  A (triangle).  In  all  cases  the  GMR  decreases  as  we  laterally 
confine  the  multilayers.  The  origin  of  this  decrease  in  the  GMR  can  be  understood 
in  terms  of  changing  the  effective  mean  free  paths  in  the  wires.  As  the  wire 
width  is  reduced  the  effective  mean  free  path  within  each  wire  decreases.  To  make 
this  more  quantitative  we  obtain  an  effective  mean  free  path  for  each  wire  and  use 
these  mean  free  paths  in  a multilayer  calculation  (infinite  width  wire).  The  results, 
which  are  shown  as  the  solid  lines,  are  in  good  agreement  with  the  exact  calculation 
(symbols).  For  cases  (a)  and  (c)  the  mean  free  paths  are  chosen  for  a Co/Cu/Co 
structure,  [181]  which  has  = 55  A,  = 10  A,  and  = 226  A.  For  case 
(b),  in  which  only  surface  scattering  contributes  to  the  GMR,  we  use  L'^  = = 

17  A and  = 226  A.  For  cases  (b)  and  (c)  the  transmission  coefficients  at  the 
interfaces  are  taken  to  be  = 0.8,  = 0.4,  [182]  while  for  case  (a),  where  the 

GMR  is  due  only  to  bulk  scattering,  we  take  = 0.6.  The  thicknesses  for 

wire  2 are  chosen  to  be  that  of  Cu  when  the  Co  slabs  are  antiferromagnetically 
coupled,  [9]  and  for  simplicity  we  choose  the  Co  layers  to  have  the  same  thickness 
as  the  Cu.  In  all  cases  the  sides  and  the  interfaces  have  a specularity  factor  S = 
0.9. 
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vs.  width  curves.  We  have  also  tried  the  same  rescaling  with  other  values  of  S, 
which  is  still  the  same  for  all  sides  and  interfaces,  and  found  the  same  curve. 
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Figure  2.14:  Rescaled  giant  magnetoresistance  as  a function  of  width.  For  large 
width  the  GMR  for  a wire  approaches  that  for  an  ordinary  unconfined  mul- 
tilayer, GMR(film).  We  define  the  width  at  which  the  GMR  is  reduced  to 
half  of  GMR(film)  as  the  half-width.  Rescaling  the  giant  magnetoresistance  by 
GMR(film)  and  the  width  by  the  half-width,  the  GMR  vs.  width  curves  of  Fig. 
2.13  fall  onto  a single  curve.  The  symbols  are  the  same  ones  used  in  Figs.  2.13(a)- 
(c). 


In  Figs.  2.15(a)  - (c)  we  have  plotted  the  half- width  as  a function  of  the  mean 
free  paths  in  the  ferromagnet  and  the  thicknesses.  The  cases  (a)  - (c)  refer  to 
the  same  parameters  as  in  Fig.  2.13.  The  ratio  between  the  mean  free  paths 
in  the  ferromagnet  is  fixed  to  L^/L^  = 5.5  in  Figs.  2.15(a)  and  (c),  which  is 
the  same  ratio  used  in  Figs.  2.13(a)  and  (c).  As  for  Fig.  2.15(b),  the  mean 
free  paths  are  equal  in  the  surface  scattering  case  of  Fig.  2.13(b):  = L^. 
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For  bulk  scattering  anisotropies  [Fig.  2.15(a)]  the  half-width  increases  almost 
linearly  with  both  the  thickness  and  the  mean  free  path.  On  the  other  hand,  for 
surface  scattering  anisotropies  [Fig.  2.15(b)]  the  half-width  depends  primarily  on 
the  thickness  of  the  film  and  only  weakly  on  the  mean  free  paths  in  the  ferromagnet. 
When  both  surface  and  bulk  scattering  anisotropies  are  present  [Fig.  2.15(c)], 
the  dependence  of  the  half-width  on  the  mean  free  path  and  thickness  can  be 
complicated.  In  the  region  where  the  half-width  has  a peak  in  Fig.  2.15(c),  the 
GMR  has  a local  minimum  as  a function  of  the  mean  free  path.  The  origin  of 
this  local  minimum  comes  from  the  near  cancellation  of  the  GMR  from  the  bulk 
and  surface  contributions.  The  location  of  this  minimum  depends  primarily  on  the 
transmission  coefficients. 

In  general  GMR  will  decrease  when  the  multilayers  are  laterally  confined,  but 
there  are  very  special  cases  where  the  GMR  actually  increases  with  lateral  con- 
finement. The  two  ways  we  have  found  to  do  this  are  (i)  by  introducing  additional 
spin  dependence  in  the  scattering  at  the  sides  of  the  wires  and  (ii)  have  the  sides 
of  the  wires  selectively  decrease  the  resistivity  in  the  ferromagnet  relative  to  the 
paramagnet.  This  again  can  be  understood  in  terms  of  the  decrease  in  the  effective 
mean  free  paths  within  each  layer  as  we  confine  the  multilayers.  By  changing  the 
mean  free  paths  by  different  amounts,  one  can  tune  the  GMR.  As  an  example  in 
Fig.  2.16  we  have  plotted  the  GMR  vs.  width  for  two  cases  which  differ  only  by 
the  specularity  factor  at  the  sides  of  the  FM.  For  the  solid  curve  the  specularity 
factors  are  0.9  for  both  spin-up  and  spin-down  electrons,  while  for  the  dashed  curve 
the  specularity  factor  at  the  sides  for  spin-up  electrons  is  0.9  and  for  spin-down 
electrons  is  0.5.  The  GMR  actually  increases  as  one  decreases  the  width  of  the 
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Figure  2.15:  Half-width  as  a function  of  the  thickness  and  mean  free  path  in  the 
ferromagnet  {L^)  for  (a)  bulk  scattering,  (b)  surface  scattering,  and  (c)  bulk  and 
surface  scattering.  In  Fig.  2.15(a)  and  (c)  the  ferromagnetic  spin-up  and  spin-down 
mean  free  paths  are  kept  at  a fixed  ratio,  /L^  = 5.5,  which  is  the  same  ratio  used 
in  Figs.  2.13(a)  and  (c).  In  the  surface  scattering  case  of  Fig.  2.15(b)  equals 
L^.  The  half- width  depends  both  on  the  thickness  and  the  mean  free  path  for  all 
three  cases.  In  case  (a),  bulk  scattering,  the  half- width  increases  roughly  linearly 
with  the  thickness  and  the  mean  free  path,  while  in  case  (b),  surface  scattering, 
the  half-width  depends  primarily  on  the  thickness  of  the  layers.  The  general  case 
when  both  bulk  and  surface  scattering  are  important  [case  (c)]  can  lead  to  complex 
dependence  on  the  mean  free  path  and  film  thickness.  The  feature  at  small  in 
(c)  is  associated  with  the  fact  that  the  GMR  has  a local  minimum  in  this  region 
(see  discussion  in  text).  Except  for  t and  L^,  which  vary,  the  parameters  used  in 
(a),  (b),  and  (c)  are  the  same  as  in  Fig.  2.13. 
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Figure  2.16:  Giant  magnetoresistance  as  a function  of  width  for  (a)  spin- 

independent  scattering  at  the  sides  and  interfaces,  (b)  spin-dependent  scattering 
at  the  sides  of  the  ferromagnet.  For  (a),  the  GMR  decreases  with  decreasing  width 
(solid  line),  whereas  for  (b),  the  GMR  increases  with  decreasing  width  (dashed 
line).  The  increase  of  the  GMR  with  reduced  width  is  due  to  the  spin-down  mean 
free  path  decreasing  faster  than  the  spin-up  one.  For  (a),  we  choose  the  same 
parameters  as  of  Fig.  2.13,  except  the  mean  free  paths  for  the  spin-up  and  down 
electrons  in  the  ferromagnet  are  200  A and  100  A,  respectively.  For  (b),  we  use  the 
same  parameters  except  that  the  specularity  factors  for  the  sides  of  the  ferromag- 
net are  spin-dependent:  = 0.9  and  = 0.5.  Note  that  the  GMR  eventually 

goes  to  zero  for  small  enough  widths  because  the  effective  mean  free  path  in  the 
paramagnet  becomes  much  smaller  than  the  thickness  of  the  paramagnetic  layer. 


75 


sample  because  the  mean  free  path  for  spin-down  electrons  decreases  more  rapidly 
than  the  mean  free  path  for  spin-up  electrons.  Such  spin  dependent  specularity 
factors  may  occur  naturally  or  be  attainable  by  coating  the  sides  of  the  multi- 
layer. Allowing  spin-dependent  specularity  factors  at  non-transmitting  interfaces 
opens  up  the  possibility  that  one  can  create  a GMR  device  without  ferromagnetic 
conductors,  but  only  insulating  ferromagnets  which  change  the  surface  scattering. 

2.5  Conclusion 

We  have  studied  the  effect  of  lateral  confinement  on  the  giant  magnetoresistance 
in  the  CIP  geometry.  With  spin  independent  specularity  factors  at  the  sides  of 
the  wires  the  GMR  in  general  decreases  as  we  reduce  the  width  of  the  wires.  In 
this  case  we  found  that  the  GMR  vs.  width  curves  could  be  very  nearly  collapsed 
onto  a single  curve  by  rescaling  the  GMR  by  its  infinite  width  value  and  the  width 
by  the  half-width.  The  half-width  depends  on  both  the  mean  free  paths  and  the 
thickness  of  the  films.  For  the  case  of  the  GMR  due  solely  to  surface  scattering, 
the  thickness  of  the  sample  plays  a dominant  role  in  determining  the  half- width, 
while  in  the  case  of  GMR  due  solely  to  bulk  scattering,  the  half-width  increases 
roughly  linearly  with  both  the  mean  free  paths  in  the  ferromagnetic  wires  and 
the  thicknesses.  The  general  case  when  both  surface  and  bulk  anisotropies  are 
important  can  lead  to  more  complex  dependences. 

The  source  of  the  decrease  in  the  GMR  in  the  CIP  geometry  is  a reduction 
in  the  effective  mean  free  path  in  the  layers  due  to  scattering  off  the  sides  of  the 
wires.  We  showed  this  quantitatively  by  determining  an  effective  mean  free  path 
within  each  wire  and  substituting  it  into  a multilayer  calculation  for  films  (infinite 
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width  wires).  The  results  of  the  approximate  solution  and  the  exact  solution  agree 
quite  well. 

For  the  case  of  spin-dependent  mean  free  paths  we  find  that  there  are  parameter 
regimes  where  the  GMR  in  CIP  geometry  can  increase  as  the  sample  is  laterally 
confined.  The  origin  of  the  increase  is  again  changing  of  the  effective  mean  free 
path  within  each  layer  as  one  decreases  the  width.  With  spin-dependent  specularity 
factors  one  can  change  the  ratio  of  the  mean  free  paths  for  spin-up  and  spin-down 
electrons  and  hence  change  the  GMR.  Thus,  with  appropriately  prepared  sides  of 
the  wires,  one  may  be  able  to  increase  or  at  least  stem  the  decrease  in  the  GMR.  In 
any  case  the  length  scale  at  which  the  reduction  in  the  GMR  takes  place  is  typically 
quite  small,  of  order  the  mean  free  paths  and  the  thicknesses  of  the  layers.  The 
effects  discussed  here  will  not  be  important  until  one  goes  to  very  small  samples. 


CHAPTER  3 

JUNCTION  MAGNETORESISTANCE  (JMR)  IN  DOUBLE  TUNNEL 
JUNCTIONS  WITH  MAGNETIC  METALS 

Magnetoresistance  in  ferromagnetic  tunnel  junctions  have  been  studied  exten- 
sively both  experimentally  and  theoretically  in  the  last  few  years  [63,  190].  More 
recently  there  are  experiments  where  magnetoresistance  has  been  studied  in  the 
Coulomb  blockade  regime.  Experimental  results  on  these  recent  experiments  were 
shown  in  Section  1.2  of  the  introduction.  The  purpose  of  this  chapter  is  to  theoret- 
ically investigate  in  detail  the  magnetoresistance  of  a double-tunnel-junction  with 
magnetic  metals  in  the  Coulomb  blockade  regime,  i.e.  at  very  low  temperature 
where  the  charging  energy  exceeds  the  thermal  energy.  The  layout  of  the  chapter 
is  the  following.  We  first  begin  with  the  theoretical  formalism  of  the  Coulomb 
blockade.  We  recover  the  well  known  results:  the  Coulomb  gap  and  the  Coulomb 
staircase  structure  in  the  current-voltage  characteristic  curves.  Dependencies  on 
different  parameters  are  also  studied.  In  the  next  section  we  use  the  same  model 
to  calculate  the  non-linear  conductances  and  the  junction  magnetoresistance  of  a 
double-tunnel-junction  with  magnetic  metals.  The  results  are  studied  in  detail. 
Finally  in  the  last  section  we  summarize  our  main  results. 

3.1  Coulomb  Blockade  in  Double  Tunnel  Junctions 

Coulomb  blockade  in  the  current- voltage  characteristic  curves  is  a manifestation 
of  single-electron  tunneling  through  nano-scale  tunnel  junctions.  The  experimental 
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realization  of  a double-tunnel-junction  consists  of  three  metals  separated  by  oxide 
layers  which  act  as  insulators.  The  two  end  leads,  denoted  by  left  (L)  and  right 
{R),  are  connected  to  a constant  voltage  source  V.  This  double  junction  can  be 
modeled  as  a closed  loop  electrical  circuit  with  two  leaky  capacitors  of  capacitances 
Cl  and  Cr  connected  in  series  (Fig.  3.1).  The  resistances  of  the  junctions,  Rl 
and  Rr,  are  connected  in  parallel  to  Cl  and  Cr  respectively.  Depending  on  the 


Figure  3.1:  Schematic  diagram  of  a double  tunnel  junction  system.  The  left  and 
right  capacitors  (Cr  and  Cr)  are  connected  in  series  with  a voltage  source.  The 
resistances  of  the  left  and  right  junctions  are  Rr  and  Rr,  respectively.  The  voltage 
in  the  left  and  right  leads  are  kept  fixed  at  V and  0,  whereas  in  the  central  region 
the  voltage  Vm  changes. 


number  of  excess  electrons,  n,  in  the  center  metal,  the  voltage  in  the  central  region 
Vm  fluctuates.  The  voltage  drops,  [Vr  - Vw(n)]  and  [Vm(u)  - Vr],  across  the  left 
and  the  right  junctions  are 

Vr  — VM(n)  = ~C 

= ^V-^-Vp,  (3.2) 

where  — e is  the  charge  of  the  electron  and  C is  the  capacitance  (Cr  + Cr).  The 
offset  voltage,  Vp,  accounts  for  the  fact  that  even  when  n and  V are  zero,  Vm  may 
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not  be  zero  due  to  the  alignment  of  the  Fermi  energy  in  the  middle  region.  Such 
effects  are  seen  experimentally  [97,  191]  in  a shifting  of  the  I - V characteristic 
away  from  a symmetric  curve  [I{V)  — I{-V)].  A nonzero  Vp  can  also  be  created 
deliberately  by  patterning  a gate  capacitor  to  the  central  region  [192,  193,  194],  To 
describe  the  transport  of  electrons  we  use  a semiclassical  model  [83,  195,  196,  197]. 
The  different  tunneling  rates  that  enter  into  the  dynamics  are:  the  rate  for  electrons 
to  tunnel  onto  the  center  metal  from  the  left  and  right  (F^_^„_^i),  and 

the  rate  for  the  electrons  to  leave  the  center  metal  to  the  left  (F^_^„_^)  and  right 
They  are  pictorially  shown  in  Fig.  3.2.  These  tunnel  rates  are  calculated 


Figure  3.2:  The  four  different  tunneling  rates  (F’s)  are  shown  pictorially.  The 
probability  that  there  are  n electrons  at  the  center  metal  at  any  time  t is  Pn{t). 

using  Fermi’s  Golden  rule.  As  an  example,  is  calculated  by  integrating  over 

energy  {dE),  the  square  of  the  tunneling  matrix  element  coupling  the  initial  and 
final  states  [|T(F')[^],  the  number  of  occupied  initial  states,  and  the  number  of 
unoccupied  final  states: 

27T 

= y J_^dE\TiE)\'‘N^(E-EL)NM{E-EM)f(E-Ei,)ll-f{E-EM)], 

(3.3) 
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where  f{E)  = 1/[1  + exp(^[^  — £^f])]  is  the  Fermi-distribution  function,  ^ is 
related  to  the  temperature  T via  ^5  = l/k^T.  Nl{E),  Nm{E)  are  the  density  of 
states  of  the  left  and  the  center  metals,  and  E^,  Em  are  their  Fermi  energies.  For 
simplicity,  we  assume  Nl{E),  Nm(E)  and  \T{E)\‘^  are  energy  independent,  so  that 
Ni{E)  w Nl,  Nm{E)  « N^,  and  \T{E)\'^  « l^oP-  With  these  assumptions  Eq. 
(3.3)  becomes 


^ n->n+l 


StT  /*oo 

= y / dEf{E  - Bi)[l  - f{E  - Em)].  (3.4) 


We  define 


Rl  = 


h 


(3.5) 


27re=JV£JVt|T„p 

to  be  the  resistance  of  the  left  lead.  The  resistance  of  the  right  lead  can  be  defined 
in  a similar  way  by  replacing  L R.  One  of  the  basic  requirement  for  single 
charge  tunneling  phenomenon  to  occur  is  that  the  tunneling  resistances  R^,  R^ 
must  exceed  the  resistance  quantum  Rk  = h/2Tre^  ~ 25.8  kf)  i.e.,  Rl(^r)  » R^. 
Using  the  identity 


/(B)li-/(^+x)i  = ^M__|g^, 


(3.6) 


the  integral  in  Eq.  (3.4)  can  be  evaluated  analytically.  The  end  result  for  the 
tunnel  rate  is 


[1  - exp(-;5[Ei  - Em]) 


Er  - E. 


M 


(3.7) 
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Similarly  it  follows  that, 


^ n->n+l  ~ 


( ^ 1 

1 — — 
1 

\emRj 

1-&^v{-^[Er- Em])_ 

For  convenience  we  introduce  a function  7(e): 


(3.8) 


7(e) 


e 

1 — 


In  terms  of  the  above  function  the  various  tunnel  rates  are: 


(3.9) 


n->ndbl 


(3.10) 


We  now  need  to  calculate  the  difference  between  the  right  (and  left)  electrode  Fermi 
energy  before  a tunneling  event  and  the  energy  of  the  center  electrode  after  the 
event,  i.e.  (Er  - Em)  [and  {El  - Em)]  is  the  energy  the  electron  gains  during  the 
tunneling  process.  For  our  calculation  we  consider  the  voltage  source  (the  battery) 


Figure  3.3:  Schematic  diagram  showing  charge  transfer  in  a double-tunnel-junction 
to  be  a capacitor  with  a capacitance  C much  larger  than  the  capacitances  of  the 
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left  [Cl)  and  the  right  {Cr)  junctions.  In  equilibrium  we  have  from  KirchoflF’s 
voltage  law 

O Or  Or, 

(3.11) 


Q _ Qi^  Or 
C~  Cr' 


where  Ql,Qr,  and  Q are  the  charges  on  the  left,  right,  and  the  big  capacitors 
respectively.  Now  let  us  consider  that  an  electron  crosses  the  left  junction  [See 
Fig.  3.3],  changing  the  charge  in  the  central  region  from  ne  to  (n  - l)e  where  n 
is  the  excess  electrons  in  the  central  region.  The  charge  in  the  left  junction  also 
decreases  from  Qi  to  (Qi  — e)  and  the  charge  in  the  second  junction  is  constant. 
The  new  charges  [Qr  — e)  and  Qr  do  no  longer  satisfy  electrostatic  equilibrium. 
To  reestablish  equilibrium  the  big  capacitor  transfers  a charge  q to  the  left  and  the 
right  capacitors.  Thus  we  have 


+ — e + Qr  + q 


Cr 


'R 


(3.12) 


Eqs.  (3.11)  and  (3.12)  can  be  solved  for  q with  the  result 


q = 


ea 


R 


{Cl  + Cr)-ClCrIC 


(3.13) 


The  total  increase  in  energy  (AE  = Em  — Er)  due  to  this  tunneling  process  is: 


AE  = 


[Ql  - e + q]^  [QR  + q]^  [Q  + q]^ 

2Cr  2Cr  2C 

2Cr  2Cr  2C' 


(3.14) 
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Substituting  the  value  of  q from  Eq.  (3.13)  and  taking  the  limit  C » Cr,Cl  we 
obtain 


Em  — Ei,  = 


Cl  + Cr  [2 


e 

- — ne 


eVC 


R 


Cl  + Cr’ 


- - ne  - CrV 


Cl  + Cr  l 
Ec-e  [Vl  - Vm{ti)]  . 


(3.15) 


Similarly  we  can  calculate  the  change  in  the  energy  across  the  right  junction. 
Considering  both  of  them  we  have 

El-  Em  = -e[VM{n)  - Vl]  - Ec,  (3.16) 

Er-  Em  = e[VM{n)  - Vr]  - Ec-  (3.17) 


Here  Ec  — e^/2(Cx,  + Cr)  is  the  charging  energy  of  the  system.  Substituting  Eqs. 
(3.16)  and  (3.17)  in  Eq.  (3.10)  we  obtain 

The  transport  of  the  electrons  thorough  the  tunnel  junctions  is  given  by  a 
master  equation  which  describes  the  time  evolution  of  the  probability,  p„(t),  for  n 
excess  electrons  in  the  center  metal  [83] 


dPn  {i) 


dt 


[Pn+lTn+l-^n  Pnrn->n+l] 
~ [Pn^rj-^n— 1 Pn— irn-l->n] 


(3.19) 
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where  the  net  tunnel  rate  is  equal  to  Defining  a matrix  M 

with  components 


= 


if  m = n ± 1 

Tn— >n+l  rti— vn— 1 if  m = 71 

0 otherwise, 


the  master  equation  may  be  written  in  a compact  form 


Mi) 

dt 


= Mp{t). 


(3.20) 


(3.21) 


Due  to  Em-^mn  = 0 which  follows  from  Eq.  (3.20)  and  Eq.  (3.21),  p„(t)  is 
independent  of  time.  The  steady  state  solution  (dp„(t)/dt  = 0)  then  has  to  satisfy 

- 0.  (3.22) 


At  steady  state,  the  current  in  the  left  junction  equals  the  current  in  the  right 
junction  due  to  current  conservation,  and  is  given  by 


7* — UW 


(3.23) 


To  write  Eq.  (3.23)  succinctly  we  define  a new  matrix  where 


V"  = 


if  m = n - 1 
if  m = n + 1 
0 otherwise. 


(3.24) 
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The  tunneling  current  can  also  now  be  expressed  in  a compact  form  as 

m 

= -eTr{?;^p(t)}.  (3.25) 

To  evaluate  the  current  in  general,  Eqs.  (3.21)  - (3.25)  have  to  be  solved 
numerically.  We  follow  the  procedure  described  in  Ref.  [196].  We  diagonalize  the 
matrix  M and  determine  its  eigenvalues  and  eigenvectors.  The  eigenvector  with 
the  zero  eigenvalue  is  the  steady  state  solution  p°.  At  low  temperature  and  voltage, 
the  higher  charge  states  |n|  >>  1 are  energetically  forbidden  and  therefore  we  can 
neglect  them.  We  found  that  for  our  choice  of  parameters  it  is  sufficient  to  keep 
the  19  states  around  n = 0.  The  highest  and  the  lowest  states  are  n = 9 and 
n = —9. 

We  restrict  ourselves  to  low  temperature  and  in  the  regime  where  the  charging 
energy  is  larger  than  the  thermal  energy,  i.e.  ksT  « Ec-  The  number  of  inde- 
pendent parameters  in  our  model  are  the  ratio  of  the  resistances  Rl/{Rl  + Rr), 
the  ratio  of  the  capacitances  Cl/ {Cl  + Cr),  and  the  additional  voltage  Vp.  In 
Fig.  3.4  we  plot  the  tunneling  current  as  a function  of  voltage  for  CVp  = 0 and 
CVp  = 0.25  for  ksT  = O.OlEc-  As  expected  for  an  asymmetric  junctions  with 
Rl  = 0.01{Rl  + Rr)  and  Cl  = 9.91{Cl-\-Cr),  the  curves  exhibit  a step-like  struc- 
ture. This  is  known  as  Coulomb  staircase.  Also  note  that  there  is  no  current  if  the 
external  voltage  is  smaller  than  e/2C . This  shift  in  the  current-voltage  character- 
istic is  called  the  Coulomb  gap  and  the  phenomenon  of  suppression  of  the  current 
below  e/2C  is  referred  to  as  Coulomb  blockade.  Comparing  the  two  curves  in  Fig. 
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3.4,  one  with  a finite  offset  voltage  and  the  other  without  the  oflFset  voltage  we  find 
that  the  effect  of  a finite  offset  voltage  Vp  is  to  shift  the  / - ^ curves.  Hence,  we 
will  consider  = 0 for  our  later  calculations. 

To  study  the  dependence  of  the  tunneling  current  on  the  ratio  of  the  capaci- 
tances and  resistances  we  plot  in  Fig.  3.4  the  I — V curves  for  different  values  of  the 
ratios.  For  asymmetric  junctions  with  Rl  « Rr  and  Cl  « Cr  (or  Rl  » Rr 
and  Cl  » Cr)  the  Coulomb  staircase  structure  is  seen.  With  increase  in  the 
ratios  Rl/{Rl  + Rr)  and  Cl! {Cl  + Cq)  the  step-like  structure  diminishes  and 
finally  for  Rl  « Rr  and  Cl  »Cr  (or  Rl  » Rr  and  Cl  « Cr)  no  step-like 
structure  is  seen. 

The  effect  of  temperature  is  to  smear  out  the  step-like  structure  seen  in  the 
I-V  curves.  This  can  be  seen  by  comparing  Fig.  3.5  with  Fig.  3.4.  With  increase 
in  temperature  the  number  of  accessible  states  increases  which  thus  destroys  the 
step-like  pattern. 

3.2  Magnetoresistance  in  the  Coulomb  Blockade  Regime 

In  this  section  we  will  study  the  different  transport  properties  of  a double-tunnel 
junction  with  magnetic  metals,  specifically  the  conductances,  the  differential  con- 
ductances for  parallel  and  antiparallel  magnetization  alignment  and  the  magne- 
toresistance. To  estimate  the  junction  resistances  we  will  use  Julliere’s  model  [65] 
which  has  been  described  in  the  section  1.2  of  Chapter  1.  In  this  model,  for  ex- 
ample the  parallel,  and  antiparallel  conductances  for  a FM/I/FM  tunnel 


87 


CpV/e 


Figure  3.4:  Current- Voltage  curves  at  low  temperature,  ksT  = O.OlE^c  for  two 
different  values  of  the  oflFset  voltage  Vp  = 0.25  and  1^  = 0.  In  both  cases  we  see 
the  Coulomb  staircase  structure  and  the  Coulomb  gap.  The  suppression  of  the 
current  at  low  voltages  is  known  as  the  Coulomb  blockade.  As  an  example,  note 
that  for  = 0 current  is  suppressed  below  e/2C.  The  width  of  the  Coulomb  gap  is 
e/C.  A finite  Vp  only  shifts  the  I — V curves.  In  both  cases  the  ratio  of  resistances 
and  capacitances  are  kept  fixed  to  Rl/{Rl  + Rr)  = 0.01  and  Cl/{Cl  + Cr)  = 0.01. 
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CrV/g  CrV/g 


Figure  3.5:  Tunneling  current  as  a function  of  voltage  for  different  values  of  the 
ratio  of  the  resistances  and  capacitances.  For  case  (a)  we  keep  the  ratio  of  the 
capacitances  fixed  to  Ci/{Cl  + Cr)  = 0.01  and  vary  the  ratio  of  Rl/{Rl  + Rr). 
For  case  (b)  we  keep  Rl/{Rl  + Rr)  to  be  fixed  to  0.01  and  vary  the  ratio  of 
the  capacitances  Cl/{Cl  + Cr).  In  both  cases  with  increase  in  the  ratio  of  the 
resistances  or  the  capacitances  the  staircase  structure  disappears.  The  temperature 
is  kept  fixed  {ksT  = O.lEc). 
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junction  are  given  by: 


Here  the  tunneling  matrix  element  T,  is  taken  to  be  both  energy  and  momentum 
independent.  The  density  of  states  of  the  electrons  at  the  Fermi  level  in  the  two 
ferromagnets  are  Ni  and  N2  and  the  fraction  of  the  number  of  spins  parallel  to  the 
magnetization  in  the  two  FMs  are  Oi  and  02  respectively.  Using  the  definitions  of 
the  spin  polarizations  Pi  and  P2  of  the  two  FMs 


Eqs.  (3.26)  and  (3.27)  can  be  rewritten  in  terms  of  Pi  and  P2. 

For  magnetic  materials,  there  are  two  types  of  carriers:  the  majority  electrons 
with  spins  parallel  to  the  magnetization  and  the  minority  electrons  with  spins  an- 
tiparallel to  the  magnetization.  We  make  the  assumption  that  the  rate  of  tunneling 
through  the  junctions  is  much  smaller  than  the  rate  of  spin  relaxation.  Under  this 
assumption  the  number  of  majority  and  minority  electrons  in  the  center  metal  is 
not  conserved  separately,  but  their  sum  is  conserved.  The  time  evolution  of  the 
probability,  p„(t),  for  n excess  electrons  in  the  center  metal  [83]  is  described  by 
the  master  equation  [Eq.  3.21].  Following  the  same  procedure  as  described  in 
the  context  of  Coulomb  blockade  the  current  is  calculated  using  Eq.  (3.25).  The 


Pi  — 2oi  — 1, 


(3.28) 

(3.29) 


P2  — 2u2  — !• 


90 


total  conductance  G is  then  simply  given  by  dividing  the  current  with  the  voltage, 
G = I /V , and  the  differential  conductance  Gdiff  is  obtained  by  differentiating  the 
current  with  respect  to  the  voltage,  G^  = dl /dV. 

Within  the  approximations  we  made  there  are  two  possible  choice  of  materials 
which  lead  to  non-zero  magnetoresistance  in  a two  junction  system:  (1)  all  the 
leads  are  ferromagnets,  i.e.  FM/I/FM/I/FM,  and  (2)  the  left  or  the  right  lead  is 
a paramagnet,  i.e.  either  FM/I/FM/I/PM  or  PM/I/FM/I/FM.  Note  that  with 
the  center  lead  non-magnetic  (FM/I/PM/I/FM)  there  is  no  JMR  or  DJMR  be- 
cause the  antiparallel  conductance  is  identical  to  the  parallel  conductance  within 
our  assumption.  For  each  of  these  cases,  we  calculate  the  conductances  for  two 
different  relative  orientations  of  the  ferromagnet  magnetizations:  parallel  and  an- 
tiparallel. Here  parallel  (P)  means  that  the  magnetizations  of  all  ferromagnets  are 
in  the  same  direction,  and  antiparallel  (AP)  means  the  magnetization  of  the  center 
ferromagnetic  metal  has  been  reversed.  In  case  (1)  of  three  ferromagnetic  metals, 
other  configurations  are  also  possible:  the  magnetization  of  the  right  or  the  left 
lead  can  be  reversed  keeping  the  magnetizations  of  the  other  two  parallel.  For  case 
(2)  where  one  of  the  left  or  the  right  lead  is  non-magnetic,  the  magnetization  of 
right  or  the  left  ferromagnetic  lead  can  also  be  reversed.  However,  the  qualitative 
features  of  the  current,  conductance,  and  the  magnetoresistance  curves  remain  the 
same.  The  junction  magnetoresistance  (JMR)  is  defined  as  the  ratio  of  the  change 
in  the  conductance  from  parallel,  G^  to  antiparallel  alignments,  G^^,  divided  by 
the  parallel  conductance 


JMR  = 1 - 


gap 


(3.30) 
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The  differential  junction  magnetoresistance  (DJMR)  is  defined  in  the  same  fashion 
as  the  JMR  but  with  the  total  conductances  replaced  by  the  differential  conduc- 
tances, Gdiff, 

^AP 

DJMR=1--^.  (3.31) 

Gdiff 

We  now  consider  the  different  parameters  of  our  model.  They  are  the  capaci- 
tances Cl  and  Cr,  the  parallel  resistances  and  R^,  the  antiparallel  resistances 
Rl^  and  R^^,  the  temperature  T,  and  the  offset  voltage  Vp.  To  estimate  the 
antiparallel  resistances,  we  use  Julliere’s  model  for  FM/I/FM  tunnel  junctions, 
where  the  junction  magnetoresistance  is  expressed  as  a product  of  the  magnitudes 
of  the  polarizations  of  the  ferromagnets:  [65] 

AR  2P1P2 
~ (l  + PiPa)' 

Here  AR  is  the  change  in  resistance  from  antiparallel  to  parallel  orientations, 
R^^  is  the  junction  resistance  when  the  magnetizations  of  the  ferromagnets  are 
antiparallel,  and  Pi,  P2  are  the  spin  polarizations  of  the  two  ferromagnets.  Eq. 
(3.32)  can  be  obtained  from  Eqs.  (1.11)  [See  Section  1.2  of  Chapter  1 for  details]. 

In  Fig.  3.6  the  currents  and  the  total  conductances  are  plotted  for  parallel 
and  antiparallel  orientations  of  the  magnetizations.  Asymmetric  junctions  and 
low  temperature  are  used  to  obtain  the  steps  in  the  current.  As  expected  in  the 
Coulomb  blockade  regime,  the  current  increases  by  steps  (cases  a and  b),  and 
the  conductance  shows  oscillations  with  decreasing  amplitude  (case  c).  Case  (a) 
corresponds  to  a finite  offset  voltage  whereas  cases  (b)  and  (c)  correspond  to  zero 
offset  voltage.  The  effect  of  a finite  Vp  is  to  shift  the  I — V curves.  For  convenience 
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we  thus  take  the  offset  voltage  to  be  zero  for  our  later  calculations.  Finally,  for 
large  enough  voltage,  the  current  becomes  linear  with  voltage  and  the  conductance 
approaches  the  classical  limiting  value  which  is  the  inverse  of  the  resistance  of  two 
resistors  connected  in  series 


1 

{Rl  + Rr) 


(3.33) 


It  is  also  interesting  to  look  at  the  diflFerential  conductances  {dl/dV)  for  parallel 
and  antiparallel  alignments  of  magnetizations.  They  are  plotted  as  a function  of 
voltage  in  Fig.  3.7.  We  see  sharp  spikes  in  the  differential  conductance  curves  at 
places  where  the  voltage  is  equal  to  ne/2G  (n  = 1,2,3,...).  At  these  places  the 
current  also  increases  by  a step  as  seen  in  Fig.  3.6. 

The  effect  of  temperature  dependence  is  shown  in  Fig.  3.8.  With  the  rise  in 
temperature,  the  number  of  accessible  states  increases.  The  steps  in  the  current 
round  off  (case  a),  and  the  amplitude  of  the  oscillations  in  the  conductance  vs. 
voltage  curve  (case  b)  gradually  decreases.  The  peaks  in  the  JMR  (case  c)  and 
the  DJMR  vs.  voltage  curves  (case  d)  become  broader  and  reduced  in  size.  The 
reduction  of  JMR  with  increase  in  temperature  has  been  observed  recently  [198]. 

Considering  the  parameter  regime  where  ksT  « Ec  and  zero  offset  voltage 
[CVp  = 0),  in  Fig.  3.9  we  plot  the  junction  magnetoresistance  for  four  different 
cases;  (a)  FM1/I/FM1/I/FM2,  (b)  FM2/I/FM1/I/FM1,  (c)  FMi/I/FMi/I/PM, 
and  (d)  PM/I/FMi/I/FMi,  where  FMi  and  FM2  are  different  ferromagnetic  met- 
als. In  the  high  bias  regime  (the  Ohmic  regime)  for  all  four  cases,  the  current  is 
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Figure  3.6:  Current  (I)  and  total  conductance  (G)  as  a function  of  voltage,  for 
a FM/I/FM/I/FM  tunnel  junction.  For  cases  (a),  (b),  and  (c),  the  solid  lines 
represent  the  current  and  the  conductance  for  parallel  alignment  of  the  magne- 
tizations, whereas  the  dashed  lines  are  for  antiparallel  alignment.  As  a function 
of  voltage  the  current  shows  steps  (cases  a and  b),  and  the  total  conductance  os- 
cillates with  decreasing  amplitudes  (case  c).  Case  (a)  corresponds  to  finite  oflFset 
voltage  CVp  = 0.25.  Cases  (b)  and  (c)  are  for  zero  offset  voltage.  The  effect  of  a 
finite  offset  voltage  is  to  shift  the  position  of  the  / - F curves.  We  choose  the  left 
and  the  center  metals  to  be  made  of  iron  with  polarization,  Ppe  = 0.40  and  the 
right  metal  to  be  made  of  cobalt  with  polarization  Pco  = 0.34  [63].  The  different 
parameters  used  are;  Cl  = O.OIC,  = 0.01  (P^  + P^),  and  UpT  = O.OlPc- 
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Figure  3.7:  Differential  conductance  Gdicr  for  parallel  and  antiparallel  alignments 
of  magnetizations  as  a function  of  voltage.  Spikes  are  seen  when  the  voltage  is 
equal  to  ne/2C  where  n is  an  integer.  Spikes  occur  at  those  places  where  step-wise 
increase  in  current  occur.  Also  note  that  the  differential  conductance  for  parallel 
orientation  is  larger  than  the  differential  conductance  for  antiparallel  orientation. 
The  parameters  chosen  are  the  same  as  in  Fig.  3.6. 
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Figure  3.8:  Temperature  dependence  of  the  (a)  Current,  (b)  Conductances,  (c) 
JMR,  and  (d)  DJMR  for  a FM/I/FM/I/FM  tunnel  junction.  As  in  Fig.  3.6,  the 
solid  and  the  dashed  lines  represent  the  current  and  the  conductance  for  parallel 
and  antiparallel  alignments.  With  increase  in  temperature,  the  steps  in  the  current 
vs.  With  increase  in  temperature,  the  steps  in  the  current  vs.  voltage  curve  round 
off  (case  a),  resulting  in  reduced  amplitude  of  oscillations  in  the  total  conductance 
(case  b).  The  peaks  in  both  the  JMR  (case  c)  and  the  DJMR  (case  d)  curves 
broaden.  Also  the  peaks  in  the  JMR  and  the  DJMR  curves  decreases  with  increase 
in  temperature.  The  parameters  are  the  same  as  used  in  Fig.  3.6,  except  for  the 
temperature  which  is  ksT  = O.lF'c. 


96 


linearly  proportional  to  the  voltage  [196] 


[CV 


{Rl  + Rr)C 


- 1 


L e 


(3.34) 


It  then  follows  from  Eq.  (3.34),  that  the  plateaus  in  Figs.  3.9(a)  - 3.9(d)  are  given 
by  the  difference  in  the  total  resistances  from  antiparallel  to  parallel  alignments 
divided  by  the  total  resistance  for  the  antiparallel  alignment, 


JMR(plateau)  = 1 - (3-35) 

For  asymmetric  junctions  in  the  low  voltage,  low  temperature  regime  we  see 
spikes  in  the  JMR  vs.  voltage  curves  in  Fig.  3.9.  The  occurrence  of  the  first  spike 
can  be  understood  analytically.  We  confine  ourselves  to  only  three  states:  the 
n = 0 state  and  the  two  states  n = ±1  immediately  accessible  from  it.  At  very  low 
temperature,  the  state  which  is  most  likely  occupied  is  the  n = 0 state  for  Vp  = 0. 
The  tunnel  rates  Fq^i  and  Fi^o  do  not  enter  in  the  expression  for  the  tunneling 
current  at  zero  temperature,  but  at  finite  temperature  they  do  contribute.  In  the 
limit  of  zero  temperature,  where  we  replace  the  function  7(e)  in  Eq.  (3.9)  by 
7(e)  = e9{e),  the  steady  state  current  can  be  expressed  in  terms  of  the  tunnel  rates 
in  the  regime  (e/2)  < CrV  < {CR/CL){ef2)  as 


1 = 


e 


T'L  -pR 
^ -l->0 

(F^^_,  + r«i^o)‘ 


(3.36) 


In  the  derivation  of  Eq.  (3.36)  we  have  assumed  that  the  capacitance  of  the 
junction  is  larger  than  the  capacitance  of  the  left  junction,  i.e,  Cr  > Cl. 
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Near  the  special  point  CrV  « e/2,  where  the  current  increases  by  a step  at  zero 
temperature,  the  height  of  the  spike  is  simply  given  by  the  difference, 

JMRCC^V'  = T ^ 0)  = 1 - Jt.  (3.37) 

At  finite  temperature,  the  height  of  the  spikes  depend  on  the  temperature,  the 
capacitances,  and  the  resistances  of  the  junctions.  We  also  notice  that  spikes  go 
up  for  cases  3.9(a)  and  3.9(c)  whereas  spikes  go  down  for  cases  3.9(b)  and  3.9(d). 
Cases  3.9(a)  - 3.9(b)  or  cases  3.9(c)  - 3.9(d)  differ  only  by  the  exchange  of  the 
metals  in  the  left  and  the  right  leads.  In  all  the  above  cases  we  keep  the  resistance 
of  the  left  junction  for  the  parallel  alignment  fixed,  = 0.01  (i?^  + R\).  By  only 
exchanging  the  metals  of  the  left  and  the  right  leads,  we  go  from  spikes  which  go  up 
to  spikes  which  go  down  and  vice  versa.  The  interchange  of  the  spikes  from  up  to 
down  or  down  to  up  can  be  explained  in  the  following  way.  At  zero  temperature, 
whenever  the  JMR  obtained  from  Eq.  (3.37)  is  greater  than  the  height  of  the 
plateau  [Eq.  (3.35)],  we  see  an  up-spike  [Figs.  3.9(a)  and  3.9(c)],  otherwise  we 
see  a down-spike  [Figs.  3.9(b)  and  3.9(d)].  Within  Julliere’s  model  [65],  we  should 
mention  that  to  see  the  spikes,  the  left  or  the  right  metal  should  be  of  a different 
ferromagnetic  material.  When  all  three  metals  are  of  the  same  ferromagnetic 
materials,  the  steady  state  solution  p°(t)  is  the  same  for  parallel  and  antiparallel 
orientation  of  magnetizations  since  the  solution  depends  only  on  the  ratio  of  the 
left  and  the  right  resistances. 

Comparing  Figs.  3.9(a)  - 3.9(b)  with  Figs.  3.9(c)  - 3.9(d),  we  see  that  the 
variation  of  the  JMR  as  a function  of  voltage  is  larger  when  one  of  the  leads  is 
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non-magnetic.  Following  Eq.  (3.35)  and  Eq.  (3.37),  which  give  the  height  of  the 
plateau  and  the  height  of  the  first  spike  at  zero  temperature,  we  see  that  the  height 
of  the  plateau  for  case  (c)  is  lower  than  for  case  (a).  On  the  other  hand,  the  heights 
of  the  first  spikes  are  the  same  for  both  cases  because  the  metals  in  the  left  and 
the  center  are  same.  Thus  the  variation  of  the  JMR  for  case  (c)  is  more  than  case 
(a).  Cases  (b)  and  (d)  can  be  understood  similarly. 

We  plot  the  differential  junction  magnetoresistance  as  a function  of  the  voltage 
in  Fig.  3.10  for  the  same  cases  as  in  Fig.  3.9  and  in  the  same  parameter  regime. 
Many  of  the  features  of  the  DJMR  are  the  same  as  for  the  JMR.  The  spikes  are 
at  the  same  places.  The  height  of  the  spikes  is  also  the  same  as  in  Eq.  (3.37). 
The  DJMR  at  large  voltages  is  the  same  as  in  Eq.  (3.37).  The  DJMR  at  large 
voltages  is  obtained  from  the  classical  value  with  two  resistors  in  series.  Finally, 
the  variation  of  the  DJMR  with  the  voltage  is  larger  in  cases  3.10(c)  and  3.9(d) 
than  in  cases  3.9(a)  and  3.9(b).  There  are,  however,  some  differences.  The  DJMR 
can  be  negative,  as  shown  in  Fig.  3.9(c),  while  the  JMR  must  be  positive  from  Eq. 
(3.37)  as  long  as  and  Cr  > Cl-  Also,  comparing  Fig.  (3.9)  with  Fig. 

(3.10),  we  see  that  the  variation  of  DJMR  as  a function  of  voltage  is  larger  than 
the  variation  of  the  JMR  with  the  voltage. 

3.3  Conclusion 

We  have  studied  the  magnetoresistance  for  double  tunnel  junctions  with  mag- 
netic metals  in  the  Coulomb  blockade  regime.  For  both  the  parallel  and  antiparallel 
orientations  of  the  magnetizations,  the  steps  in  the  current  and  the  oscillations  in 
the  total  conductance  are  observed.  Also  with  increase  in  temperature  the  struc- 
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Figure  3.10:  Differential  Junction  magnetoresistance  (DJMR)  as  a function  of 
voltage  for  the  same  cases  as  in  Fig.  3.9.  We  see  spikes  in  the  DJMR  at  those 
places  where  the  current  increases  by  a step.  As  in  Fig.  3.9,  the  variation  in  the 
DJMR  is  larger  as  a function  of  voltage  for  cases  (c)  and  (d)  compared  to  cases 
(a)  and  (b).  For  case  (c)  the  DJMR  changes  sign.  Also  compared  to  Fig.  3.9,  the 
change  in  the  DJMR  is  larger  than  the  JMR.  The  metals  are  the  same  as  used  in 
Fig.  3.9  and  the  parameters  are  identical  to  Fig.  3.6. 
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ture  is  reduced.  These  are  expected  in  the  Coulomb  blockade  regime.  The  JMR  as 
a function  of  voltage  shows  spikes  at  the  same  places  where  the  current  increases 
by  a step.  Although  the  height  of  the  spikes  in  general  depends  on  the  tempera- 
ture, capacitances,  and  the  resistances  of  the  system  the  the  height  of  the  plateaus 
in  the  JMR  vs.  voltage  curves  are  given  by  the  classical  value  obtained  by  adding 
two  resistors  in  series.  For  asymmetric  junctions  where  Cr  > Cl,  the  spikes  go 
up  if  the  ratio  of  the  resistances  of  the  right  junction  for  parallel  and  antiparallel 
alignment  {R^/Rr^)  is  greater  than  the  same  ratio  of  the  resistances  of  the  left 
junction  (R^/R^^).  Otherwise,  the  spikes  go  down. 

Many  of  the  features  of  the  DJMR  are  similar  to  the  JMR:  the  spikes  occur  at 
the  same  places,  the  height  of  the  first  spike  at  zero  temperature  is  the  same,  and 
the  sign  (up  or  down)  of  the  spikes  are  the  same.  Also  at  high  bias,  the  DJMR 
saturates  to  the  classical  value  obtained  by  adding  two  resistors  in  series.  However, 
there  are  few  differences.  The  DJMR  can  be  negative  for  some  cases,  while  the 
JMR  is  positive  as  long  as  < R^^  and  Cr>  Cl-  Finally,  the  variation  of  the 
DJMR  is  larger  as  a function  of  voltage  than  that  of  the  JMR  . 

Thus  double  tunnel  junction  systems  which  have  magnetic  metals  exhibit  rich 
Coulomb  blockade  conductance  vs.  voltage  curves.  These  curves  provide  a signa- 
ture of  both  the  Coulomb  charging  effects  and  the  spin  polarization  of  the  tun- 
neling. By  measuring  the  magnetoresistance,  one  can  extract  information  about 
both  the  capacitance  charging  energies  and  the  spin  polarization  of  the  tunneling 
electrons.  Hence  theories  of  spin  polarized  tunneling  and  the  Coulomb  blockade 
both  can  be  tested.  Furthermore,  there  are  cases  where  the  magnetoresistance  can 
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be  dramatically  enhanced  near  one  of  the  Coulomb  blockade  steps,  meaning  that 
there  may  be  some  applications  of  these  eflFects  in  electronic  devices  of  the  future. 


CHAPTER  4 

MAGNETIZATION  REVERSAL  IN  SINGLE-DOMAIN  THIN  FILMS 
BY  EXTERNAL  ELECTRIC  FIELD 

Magnetization  reversal  in  thin-films  by  applying  a magnetic  field  is  a well- 
known  fact.  This  phenomenon  has  been  widely  used  in  industry.  But  the  idea  that 
it  may  be  possible  to  change  the  magnetization  by  applying  an  external  electric 
field  in  some  spin-ordered  materials  was  due  to  Landau  and  Lifshitz  [171].  This 
effect  is  named  the  magnetoelectric  (ME)  effect  and  has  been  discussed  in  Section 
1.3  of  Chapter  1.  The  ME  effect  in  single-domain  ferromagnetic  thin  films  is  rather 
new.  We  propose  that  it  is  indeed  possible  that  an  electric  field  can  change  the 
magnetic  surface  anisotropy  of  a thin  film  which  thereby  changes  the  magnetization 
direction.  The  role  of  magnetic  surface  anisotropy  also  has  been  discussed  in  the 
introduction.  To  show  the  effect  of  the  electric  field  on  the  surface  anisotropy  we 
do  a simple  atomic  calculation  which  serves  as  a motivation  and  also  sets  the  stage 
for  studying  real  systems.  In  real  systems  (thin  films  of  transition  metals)  the 
applied  electric  field  will  be  screened  within  a layer.  For  real  systems  one  also  has 
to  include  the  details  of  the  band-structure  of  the  materials  [199,  200,  201]. 

4.1  Effect  of  Electric  Field  on  Magnetic  Anisotropy  — An  Atomic 

Calculation 

We  begin  with  a simple  atomic  calculation  to  evaluate  the  magnetic  surface 
anisotropy  for  a single  atom  (of  Fe,  Co,  and  Ni).  The  Hamiltonian  consists  of  three 
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parts:  Hq  which  is  the  unperturbed  part  of  the  Hamiltonian,  i^LS  which  contains 
the  spin-orbit  coupling  term,  and  Uappi  which  is  due  to  the  external  electric  field. 

H = H(j  + Hi^s  + Ua,pp\.  (4-1) 

The  spin-orbit  part  of  the  Hamiltonian  is 


Hls  = gL-S, 


(4.2) 


where  g is  the  local  spin-orbit  coupling  constant. 

The  part  of  the  Hamiltonian  due  to  the  applied  electric  field  (along  the  z 
direction)  is: 

= \e\£z.  (4.3) 

To  construct  the  full  Hamiltonian  matrix  H the  basis  functions  are  the  or- 
thonormal atomic  orbitals  \fj,a  > where  fj,  denotes  one  of  the  atomic  orbital  and  a is 
the  spin  (=t,  i).  The  atomic  orbitals  are  the  3d  subbands  which  are  yz,  zx,  xy,  di  = 
(x^  — y2),  dj  = {3z^  — r^),  the  4s  orbital,  and  the  Ap  orbitals  Px,Py,  and  p^.  x,  y,  z 
refer  to  the  standard  axes  and  the  spin  quantization  is  along  an  axis  z'  character- 
ized by  the  standard  angular  co-ordinates  Q and  <^.  They  are  shown  in  Fig.  4.1. 


The  atomic  wave  functions  for  different  orbitals  are: 
3d  orbitals: 


d 


f{r)  sin  9 cos  9 sin  (j), 


(4.4) 
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Figure  4.1:  Pictorial  representation  of  the  angles.  Magnetization  (M)  makes  an 
angle  0 with  the  z axis.  Due  to  azimuthal  symmetry,  the  magnetization  and  the 
total  energy  do  not  depend  on  the  angle  0. 


4g  orbital: 


4p  orbitals: 


dzx  = 

/l5 

\ — /(r)  sin  6 cos  6 cos  (j>, 
V 47T 

(4.5) 

dxy  = 

/ 15 

\/ — /(r)  sin^  6 sin  0 cos  (f), 
V 47T 

(4.6) 

dx2^y2  -- 

^16^/(^)sin^^cos  2(f>, 

(4.7) 

dsz^-r^  = 

\/i6^/W(3cos2  0 1), 

(4.8) 

s = 


(4.9) 


lY , , , 

Px  = 

y —h{r)  sin  9 cos  (f), 
V 47T 

(4.10) 

Py  = 

lY 

\ -—h{r)  sin  9 sin  <j>, 
V 47T 

(4.11) 

Pz  = 

[Y 

\ h(r)  cos  9, 

V 47t 

(4.12) 
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where  f{r),g{r),  and  h{r)  are  radial  part  of  the  orbitals  and  they  are  normalized 
as 


jdrr^lf{r)f  = 1, 

ldrr^[g(r)f  = 1, 

[ drr^  [/i(r)]^  = 1. 


(4.13) 

(4.14) 

(4.15) 


Using  the  atomic  wave  functions  in  Eqs.  (4.4)  - (4.12)  the  spin-orbit  matrix 
elements  < g,a\L  ■ > can  be  obtained.  In  matrix  notation  the  Hamiltonian 

in  Eq.  (4.2)  takes  the  following  form: 


u-n 

■^LS 

•^LS 

M 

N 

^LS 

TT-U 

■“LS 

-N* 

M* 

(4.16) 


The  matrix  elements  of  M and  N are  given  in  the  Appendix  A.  Similarly  the 
matrix  elements  for  are  calculated  and  they  are  also  given  in  Appendix  B. 

We  now  diagonalize  the  18  x 18  matrix,  fill  up  the  states  using  Hund’s  rule  [202] 
and  obtain  the  total  energy  of  the  filled  states  (Etot)-  In  the  absence  of  electric  field 
(crystal  or  applied)  our  atomic  Hamiltonian  in  Eq.  (4.1)  is  rotationally  invariant 
- thus  there  is  no  preferential  direction  for  the  magnetization.  By  applying  an 
electric  field  we  partially  break  the  rotational  symmetry  (the  azimuthal  symmetry 
is  not  broken)  and  the  magnetization  then  prefers  to  be  in  a particular  direction. 
To  mention  that  in  real  systems  even  in  the  absence  of  the  applied  electric  field 
the  magnetization  has  a preferred  direction.  For  our  simple  atomic  calculation  in 
the  presence  of  the  applied  field  the  total  energy  depends  only  on  the  angle  6. 
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The  change  in  the  total  energy  of  the  system  is  calculated  by  calculating  the 
total  energy  of  the  system  for  two  cases;  one  with  the  electric  field  and  the  other 
without  the  electric  field.  The  difference  of  these  two  is  the  change  in  energy  due 
the  electric  field. 

A£'  = — E^o^(£  — 0).  (4.17) 

The  different  input  parameters  of  our  problem  are  the  ground  state  energies 
for  the  s,p  and  d atomic  orbitals  {E,,  Ep,  E^,  E^)  [203],  the  exchange  energies 
(A)  between  up  and  down  spin  d-shell  conduction  electrons  [147],  the  spin-orbit 
coupling  parameter  (g)  [164],  the  atomic  spacing  (a)  [199],  and  the  two  unknown 
parameters  {p  and  q)  which  are  the  overlaps  of  the  radial  parts  of  the  3d  with  4s 
and  4p  orbitals. 

P = J drr^f{r)g{r),  (4.18) 

q — j drr^h{r)g{r).  (4.19) 


To  get  an  estimate  of  the  parameters  p and  q we  use  the  one  electron  atom  radial 
wave-functions  [202,  204] 


9{r) 

h{r) 

/(^) 


(^) 


1 - 


Q—Zr/Zas 


4v^ 

/Zr  \ 

[-(a 

3 

yzas) 

kZub  ) 

2v^ 

„-Zrl3aB 

3^5 

KZae) 

\3cib  j 

^-ZrjZaB 


(4.20) 

(4.21) 

(4.22) 
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where  Z is  the  atomic  number  and  Ob  is  the  Bohr  radius.  Using  Eqs.  (4.20  - 4.22), 
Eqs.  (4.18)  and  (4.19)  yield  after  integration 


MD  ■ 

(4.23) 

(4.24) 

We  calculate  the  anisotropy  energy  for  Fe,  Co,  and  Ni  atoms.  The  different 
parameters  chosen  are  listed  in  Table  4.1. 

Table  4.1:  Different  Input  Parameters 


Element 

E,(eV) 

B,(eV) 

Bj(eV) 

Bi(eV) 

A(eV) 

ff(eV) 

a(A) 

Fe 

-7.95 

-5.51 

-7.02 

-5.22 

1.8 

0.075 

2.861 

Co 

-8.37 

-5.38 

-7.40 

-6.40 

1.0 

0.085 

3.281 

Ni 

-8.80 

-5.48 

-8.70 

-8.20 

0.5 

0.105 

3.523 

To  show  how  the  energy  varies  with  the  electric  field  we  plot  the  change  in  the 
energy  of  a Fe  atom  in  Fig.  4.2.  The  change  in  energy  is  negative  and  also  varies 
quadratically  with  the  field.  Similar  results  are  seen  with  Co  and  Ni  atoms.  The 
quadratic  dependence  with  the  electric  field  and  the  negative  sign  is  also  expected 
for  the  energy  from  second  order  perturbation  theory. 

The  effect  of  the  electric  field  on  the  anisotropy  energy  for  Fe  and  Co  atoms  is 
shown  in  Fig.  4.3.  For  a fixed  value  of  the  electric  field  the  change  in  the  energy 
(considering  only  the  anisotropy  part  AE)  is  plotted  as  a function  of  cosine  of  the 
angle  Q.  The  anisotropy  energy  takes  the  form: 


AE  = Ki  cos^  9 + K2  cos'^  9 + ... 


(4.25) 
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in  10"^  (V/a.uf 


Figure  4.2;  Change  in  the  total  energy  of  a Fe  atom  due  to  the  electric  field 
[AEt  — E{S,6)  — E(£  = 0,^)]  is  plotted  as  a function  of  square  of  the  field  for 
two  different  values  of  the  angle  9 = 0 and  tt/2.  The  applied  field  is  in  volts 
(V)/atomic  unit  (a.u).  The  change  in  energy  is  negative  and  scales  quadratically 
with  the  field. 


no 


The  anisotropy  constants  Ki  and  K2  depend  on  the  spin-orbit  coupling  parameter 


COS  0 


Figure  4.3:  The  anisotropy  energy  as  a function  of  cos^  6 for  Fe  and  Co  atoms.  In 
all  two  cases  the  energy  is  minimum  for  0 = 0 which  implies  that  the  electric  field 
prefers  the  magnetization  to  be  perpendicular  to  the  x — y plane.  The  electric  field 
is  fixed  to  = 0.01  V/a.u  in  our  calculation  for  all  cases. 


and  also  on  the  electric  field.  K2  is  at  least  an  order  of  magnitude  smaller  than  Ki. 
The  anisotropy  energy  is  minimum  when  ^ = 0 - thus  the  anisotropy  energy  due  to 
electric  field  and  spin-orbit  coupling  prefers  the  magnetization  to  be  perpendicular 
to  the  x — y plane  (See  Fig.  4.1). 

On  the  other  hand,  for  Ni  atom  the  anisotropy  energy  due  to  the  electric  field 
and  spin-orbit  coupling  is  minimum  for  9 — tt/2.  Magnetization  in  this  case  prefers 
to  be  in  the  x — y plane.  This  is  shown  in  Fig.  4.4. 


Anisotropy  energy  (meV) 
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Figure  4.4:  The  anisotropy  energy  as  a function  of  cos^  9 for  a Ni  atom.  In  this  case 
the  energy  is  minimum  for  0 = 7t/2.  The  electric  field  thus  prefer  the  magnetization 
to  be  in  the  x — y plane.  The  electric  field  is  fixed  to  the  same  value  as  in  Fig.  4.3 
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4,2  Conclusion 

This  simple  atomic  calculation  shows  that  the  electric  field  can  indeed  effect 
the  magnetization.  As  mentioned  in  Section  1.3  of  Chapter  1,  for  real  magnetic 
systems  (thin  films  of  transition  metals)  the  magnetization  can  be  perpendicular  to 
the  plane  (for  bcc  Fe  001)  or  in-plane  (for  Co  and  Ni  films).  For  bcc  Fe  films  larger 
than  five  monolayers  a transition  from  perpendicular  to  in-plane  has  been  seen 
experimentally  [140,  141].  Our  atomic  calculation  shows  that  it  may  be  possible  to 
control  the  transition  by  applying  an  external  electric  field.  This  new  phenomenon 
holds  the  promise  of  being  applied  in  the  future  recording  devices. 


CHAPTER  5 
CONCLUSION 


In  this  thesis  I studied  transport  properties  of  three  different  types  of  nanostruc- 
tures: laterally  confined  magnetic  multilayers  (Chapter  2),  double-tunnel-j unctions 
made  with  magnetic  metals  (Chapter  3),  and  thin  ferromagnetic  films  (Chapter 
4).  This  chapter  will  review  the  major  results  obtained  and  provide  suggestions 
for  future  theoretical  work. 

5.1  GMR  in  Laterally  Confined  Magnetic  Multilayers 

In  magnetic  multilayers  I studied  the  effect  of  lateral  confinement  on  the  giant 
magnetoresistance  (GMR)  in  the  current-in-plane  geometry.  GMR  decreases  as  we 
decrease  the  width  of  the  wire  for  spin-independent  specularity  factors  at  the  sides 
of  the  wires.  This  reduction  in  GMR  can  be  understood  in  terms  of  the  reduction 
in  the  effective  mean  free  path  in  the  layers  due  to  scattering  at  the  sides  of  the 
wires.  The  length  scale  where  the  GMR  decreases  depends  on  both  the  mean-free 
paths  and  the  thickness  of  the  films.  With  spin-dependent  specularity  factors  at 
the  sides  I find  that  there  are  parameter  regimes  where  GMR  can  increase  with 
lateral  confinement. 
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5.2  JMR  in  Tunnel  Junctions  with  Magnetic  Metals 

Tunnel  junctions  made  with  magnetic  metals  exhibit  both  the  effects  of  spin- 
polarized  tunneling  and  the  phenomenon  of  the  Coulomb  blockade.  The  purpose 
of  my  study  was  to  investigate  in  the  Coulomb  blockade  regime  the  current,  con- 
ductance, and  the  junction  magnetoresistance  (JMR)  in  double-tunnel  junctions 
made  with  magnetic  metals.  I find  that  the  current-voltage  characteristic  curves 
show  a stair-case  structure  as  expected  at  low  voltage  and  low  temperatures  (the 
Coulomb  blockade  regime),  the  conductance  curves  show  oscillations  which  even- 
tually saturate  with  increasing  voltage.  On  the  other  hand,  the  JMR  as  a function 
of  voltage  exhibits  spikes  - the  spikes  may  go  up  or  may  go  down  depending  on 
the  magnetic  material  in  the  leads.  The  spikes  occur  at  those  places  where  the 
current  increases  by  a step.  At  large  voltages  (Ohmic  region)  no  spikes  are  seen  in 
the  JMR  vs.  voltage  curves. 

Recently  the  effect  of  barrier  impurity  atoms  on  the  magnetoresistance  of  fer- 
romagnetic tunnel  junctions  has  been  investigated  [205].  A submonolayer  amounts 
of  Co,  Pd,  Cu,  or  Ni  is  incorporated  into  the  middle  of  the  insulating  oxide  of 
Co/Al203/NigoFe2o  junctions.  It  has  been  observed  that  the  JMR  decays  approx- 
imately linearly  with  increasing  impurity  content  (except  for  Co  at  300  K in  Fig. 
5.1b).  The  slope  depends  sensitively  on  the  type  of  element  but  the  decrease  of 
JMR  is  much  more  pronounced  in  case  of  Ni  impurities.  This  is  shown  in  Fig. 
5.1.  Nevertheless  the  shape  of  the  JMR  curves  remains  practically  unchanged.  For 
Ni  impurities  this  is  illustrated  in  Fig.  5.2.  Although  spin-flip  scattering  may  be 
a reason  for  such  a decrease  of  JMR  with  impurity  thickness  for  most  cases  but 
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(a) 


Thickness  (A) 


(b) 


Thickness  (A) 


Figure  5.1:  Normalized  JMR  as  a function  of  the  thickness  of  the  layer  of  impurities 
present  in  the  tunneling  barrier  [205]  at  two  different  temperatures  (a)  77  K and 
(b)  300  K.  JMR  approximately  decreases  linearly  with  the  thickness.  In  case  of 
Co  at  300  K we  note  a clear  deviation  from  a linear  fit. 
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magnetic  field  (Oe) 


Figure  5.2;  Magnetoresistance  curves  for  tunnel  junctions  for  different  layers  of 
Ni  impurities  at  77  K [205].  The  shape  of  the  JMR  curves  remains  practically 
unchanged  with  the  thickness  of  Ni  impurities. 
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it  fails  to  explain  the  results  for  Co  at  300  K.  A fully  satisfactory  theory  is  still 
lacking. 


5.3  Magnetocrystalline  Anisotropy  in  Thin  Magnetic  Films 

Magnetocrystalline  anisotropy  originates  due  to  the  coupling  between  the  spin 
and  the  angular  momentum  of  the  electrons.  At  the  surface  of  the  thin  films  the 
anisotropy  is  much  more  pronounced  than  the  bulk.  This  is  known  as  the  surface 
anisotropy  and  has  been  under  intense  investigation  recently  by  both  theorists 
and  experimentalists.  Applying  an  electric  field  in  these  thin  films  can  change  the 
surface  anisotropy  and  thus  change  the  magnetization.  In  Chapter  4 with  a simple 
atomic  calculation  I showed  the  effect  of  the  electric  field  on  the  anisotropy.  Band- 
structure  calculations  for  real  systems  are  needed  to  get  quantitative  values  for  the 
anisotropy.  For  thin  films  of  transition  metals  magnetic  moment  is  enhanced  at 
the  surface  than  the  bulk.  Thus,  one  can  also  investigate  the  relationship  between 
the  surface  magnetic  moment  with  the  surface  anisotropy. 

A new  interesting  phenomena  may  be  associated  with  thin  magnetic  films. 
Pressure-induced  quantum  phase  transition  in  heavy-fermion  metals  has  been  re- 
ported [206,  207].  Quantum  critical  behavior  of  itinerant  ferromagnets  where  a 
quantum  phase  transition  from  a ferromagnetic  to  a paramagnetic  state  occurs 
has  also  been  studied  theoretically  [208].  If  we  can  prepare  a thin  magnetic  film 
near  the  critical  regime  then  by  applying  a small  electric  field  it  may  be  possible 
to  go  from  ferromagnetic  to  paramagnetic  state.  Applying  an  electric  field  changes 
the  number  density  of  electrons  in  each  layer  of  the  film  which  may  lead  to  a quan- 
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turn  phase  transition.  Experimental  studies  as  well  as  theoretical  investigations 
may  unravel  this  interesting  phenomenon. 


APPENDIX  A 

THE  SPIN-ORBIT  COUPLING  MATRIX  ELEMENTS 


As  described  in  the  text  the  spin-orbit  part  of  the  Hamiltonian  for  a monolayer 
is  of  the  form: 


-^mono  — 9 


M N 
-N*  M* 


(A.l) 


where  M and  N are  9x9  matrices  and  g is  the  constant  spin-orbit  coupling 
parameter. 

The  matrix  elements  for  M are:  {di  = di2_j,2,  d.2  = 


< Pi  t |L  • S|p2 1>  = ^sin^sin</>,  (A.2) 

<Pit  |L-S|p5,t>  = -^cos0,  (A.3) 

< Pi  t |L  • S|pj,  t>  = ^Sin0cos0,  (A. 4) 

< di  t |L  • S|cZii  t>  = ^sin0sin0,  (A.5) 

< di  t |L  • S|dj,i  t>  = ]^sin9cos(j),  (A.6) 

< di  t |L  • S|4j,  t>  = -i  cos  e,  (A.7) 

■\/3i 

< dii  t |L  • S|d2 1>  = — — sin^sin0,  (A. 8) 

A 

dxz  T ' S|dyi  'f>  = — — cos^,  (A. 9) 

< dii  t |L  • S|diy  t>  = ^sin0cos0,  (A.  10) 

< c?2 1 |L  • S|dj,i  t>  = ^^sin^cos0,  (A. 11) 
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< dy,  t |L  • S|4y  t> 


— -sin^sin(/). 


(A.12) 


Other  non-zero  matrix  elements  of  M are  obtained  by  taking  the  complex 
conjugate: 


< /3  t |L  • S|a  t>=  [<  a t |L  • S|^  t>]*.  (A.13) 


The  matrix  elements  for  N are: 

< Px  t |L  • S|px  4>  = 

< Px  t |L  • S|pj,  4,>  = 

< Px  t |L  • s|pj,  ;>  = 

< di  t |L  • S|4x  i>  = 

< di  t |L  • S|dyx  i>  = 

< di  t |L  • S|4y  i>  = 

< 4z  t |L  • s|d2  ;>  = 
dxz  T ■ S|dyx  ~ 
dxz  T ‘ S|dxy  ~ 

< d2  T |L  • s|d,x  ;>  = 

dyz  |li  • S|dxj/  ~ 


- (cos  4>  + i sin  <j)  cos  9) , 

(A.  14) 

i 

- sin  9, 

(A.15) 

— -(sin0  — zcos^cos  (j)), 

(A.16) 

-(cos  (f)  + i cos  9 sin  ^), 

(A.17) 

— ~{s\n(t>  — icos9cos  4>), 

(A.18) 

i sin  9^ 

(A.19) 

(cos  4>  + i sin  </>  cos  9) , 

(A.20) 

i 

-s\n9, 

(A.21) 

—-{sm(f)  — i cos  ^ cos  4>), 

(A.22) 

— -^^^(sin0  — i cos  cos  9), 

(A.23) 

— -(cos(^  -1-  z cos  0 sin  <j)). 

(A.24) 

Other  non-zero  matrix  elements  of  N are  obtained  by 


< /3  t |L  • S|a  ;>=  - < a t |L  • S|^  i>  . 


(A.25) 


APPENDIX  B 

ELECTRIC  FIELD  MATRIX  ELEMENTS 
The  part  of  the  Hamiltonian  which  depends  on  the  electric  field  is: 

= \e\Sz.  (B.l) 

The  matrix  elements  are  calculated  using  the  atomic  orbitals  given  in  Eqs.  (4.4)  - 
(4.12).  The  matrix  elements  also  do  not  depend  on  spin  of  the  electrons,  i.e. 


< a \e£ z\P '[>  = 

< a 4-  \e£z\l3  |>, 

(B.2) 

< a t \e£z\l3  = 

< a 4-  \e£z\l3  t>=  0. 

(B.3) 

and  they  are: 

< s\eSz\p2  >=<  Pz\e6z\s  > 

< dyz\e£z\py  >=<  Py\e£z\dy2  > 

< dxz\e£z\px  >=<  pj.\e£z\d,;z,  > 

< d2\e£z\p;z  >—<  p2\e£z\d2  > 


(B.4) 

(B.6) 

(B.6) 

(B.7) 
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where  the  two  unknown  parameters  p and  q are  overlaps  of  the  radial  parts  of  the 
3d  with  4s  and  4p  orbitals. 


P = j drr^f{r)g(r),  (B.8) 

q = J drr^h(r)g{r).  (B.9) 
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